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Abstract 

The general features of the expansion in statistical mechanics and quantum field 
theory are briefly reviewed both from the theoretical and from the phenomenological 
point of view as an introduction to a more detailed analysis of the large- properties 
of spin and gauge models possessing the symmetry group SU(A^) x SU(A^). 

An extensive discussion of the known properties of the single-link integral (equiv- 
alent to YM2 and one-dimensional chiral models) includes finite-A^ results, the ex- 
ternal field solution, properties of the determinant, and the double scaling limit. 

Two major classes of solvable generalizations are introduced: one-dimensional 
closed chiral chains and models defined on a d— 1 dimensional simplex. In both cases 
large-A solutions are presented with emphasis on their double scaling properties. 

The available techniques and results concerning unitary-matrix models that corre- 
spond to asymptotically free quantum field theories (two-dimensional chiral models 
and four-dimensional QCD) are discussed, including strong-coupling methods, re- 
duced formulations, and the Monte Carlo approach. 
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1 Introduction 

1.1 General motivation for the 1/N expansion 

The approach to quantum field theory and statistical mechanics based on the 
identification of the large- limit and the perturbative expansion in powers 
of where is a quantity related to the number of field components, is 
by now almost thirty years old. It goes back to the original work by Stanley 
[1] on the large- iV limit of spin systems with O(A^) symmetry, soon followed 
by Wilson's suggestion that the expansion may be a valuable alternative 
in the context of renormalization-group evaluation of critical exponents, and 
by 't Hooft's extension [2] to gauge theories and, more generally, to fields 
belonging to the adjoint representation of SU(A^) groups. More recently, the 
large-iV limit of random-matrix models was put into a deep correspondence 
with the theory of random surfaces, and therefore it became relevant to the 
domain of quantum gravity. 

In order to understand why the 1/N expansion should be viewed as a funda- 
mental tool in the study of quantum and statistical field theory, it is worth 
emphasizing a number of relevant features: 

1) is an intrinsically dimensionless parameter, representing a dependence 
whose origin is basically group-theoretical, and leading to well-defined field 
representations for all integer values, hence it is not subject to any kind of 
renormalization; 

2) does not depend on any physical scale of the theory, hence we may expect 
that physical quantities should not show any critical dependence on A^ (with 
the possible exception of finite- A^ scaling effects in the double-scaling limit); 

3) the large-A^ limit is a thermodynamical limit, in which we observe the 
suppression of fluctuations in the space of internal degrees of freedom; hence we 
may expect notable simplifications in the algebraic and analytical properties 
of the model, and even explicit integrability in many instances. 

Since integrability does not necessarily imply triviality, the large- A^ solution to 
a model may be a starting point for finite- A^ computations, because it shares 
with interesting finite values of A^ many physical properties. (This is typically 
not the case for the standard free-field solution which forms the starting point 
for the usual perturbative expansions.) Moreover, for reasons which are clearly. 
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if not obviously, related to the three points above, the physical variables which 
are naturally employed to parameterize large- results and expansions 
are usually more directly related to the observables of the models than the 
fields appearing in the original local Lagrangian formulation. 

More reasons for a deep interest in the study of the large- expansion will 
emerge from the detailed discussion we shall present in the rest of this in- 
troductory section. We must however anticipate that many interesting review 
papers have been devoted to specific issues in the context of the large- limit, 
starting from Coleman's lectures [3], going through Yaffe's review on the rein- 
terpretation of the large- limit as classical mechanics [4], Migdal's review on 
loop equations [5], and Das' review on reduced models [6], down to Polyakov's 
notes [7] and to the recent large commented collection of original papers by 
Brezin and Wadia [8], not to mention Sakita's booklet [9] and Ma's contri- 
butions [10,11]. Moreover, the 1/N expansion of two-dimensional spin models 
has been reviewed by two of the present authors a few years ago [12]. As a 
consequence, we decided to devote only a bird's eye overview to the general 
issues, without pretension of offering a self-contained presentation of all the 
many conceptual and technical developments that have appeared in an enor- 
mous and ever-growing literature; we even dismissed the purpose of offering a 
complete reference list grouped by arguments, because the task appeared to 
be beyond our forces. 

We preferred to focus on a subset of all large- A^ topics, which has never been 
completely and systematically reviewed: the issue of unitary-matrix models. 
Our self-imposed limitation should not appear too restrictive, when consider- 
ing that it still involves such topics as U(A^) x U(A^) principal chiral models, 
virtually all that concerns large- A?" lattice gauge theories, and an important 
subset of random-matrix models with their double-scaling limit properties, 
related to two-dimensional conformal field theory. 

The present paper is organized on a logical basis, which will neither neces- 
sarily respect the sequence of chronological developments, nor it will keep the 
same emphasis that was devoted by the authors of the original papers to the 
discussion of the different issues. 

Sect. 2 is devoted to a presentation of the general and common properties of 
unitary-matrix models, and to an analysis of the different approaches to their 
large- A?" solution that have been discussed in the literature. 

Sect. 3 is a long and quite detailed discussion of the most elementary of all 
unitary-matrix systems. Since all essential features of unitary-matrix models 
seem to emerge already in the simplest example, we thought it worthwhile to 
make this discussion as complete and as illuminating as possible. 

Sect. 4 is an application of results obtained by studying the single-link prob- 
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lem, which exploits the equivalence of this model with lattice YM2 and prin- 
cipal chiral models in one dimension. 

Sect. 5 is devoted to a class of reasonably simple systems, whose physical 
interpretation is that of closed chiral chains as well as of gauge theories on 
polyhedra. 

Sect. 6 presents another class of integrable systems, corresponding to chiral 
models defined on a d-dimensional simplex, whose properties are relevant both 
in the discussion of the strong-coupling phase of more general unitary-matrix 
models and in the context of random-matrix models. 

Sect. 7 deals with the physically more interesting applications of unitary- 
matrix models: two-dimensional principal chiral models and four-dimensional 
lattice gauge theories, sharing the properties of asymptotic freedom and "con- 
finement" of the Lagrangian degrees of freedom. Special issues, like numerical 
results and reduced models, are considered. 



1.2 Large N as a thermodynamical limit: factorization 

As wc already mentioned briefly in the introduction, one of the peculiar fea- 
tures of the large- TV limit is the occurrence of notable simplifications, that 
become apparent at the level of the quantum equations of motion, and tend 
to increase the degree of integr ability of the systems. These simplifications 
are usually related to a significant reduction of the number of algebraically- 
independent correlation functions, which in turn is originated by the property 
of factorization. 

This property is usually stated as follows: connected Green's functions of quan- 
tities that are invariant under the full symmetry group of the system are sup- 
pressed with respect to the corresponding disconnected parts by powers of 
Hence when iV — > 00 one may replace expectation values of products of 
invariant quantities with products of expectation values. 

One must however be careful, since factorization is not a property shared by 
all invariant operators without further qualifications. In particular, experience 
shows that operators associated with very high rank representations of the 
symmetry group, when the rank is 0{N), do not possess the factorization 
property. A very precise characterization has been given by Yaffe [4], who 
showed that factorization is a property of "classical" operators, i.e., those 
operators whose coherent state matrix elements have a finite N ^ 00 limit. 

It is quite interesting to investigate the physical origin of factorization. The 
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property 



lim {AB) = {A) {B) 



(1.1) 



implies in particular that 

hm{A') = {Af, (1.2) 

i.e., the vacuum state of the model, seen as a statistical ensemble, seems to 
possess no fluctuations. To be more precise, all the field configurations that 
correspond to a nonvanishing vacuum wavefunction can be related to each 
other by a symmetry transformation. This residual infinite degeneracy of the 
vacuum configurations makes the difference between the large- limit and 
a strictly classical hmit — > 0, and allows the possibility of violations of 
factorization when infinite products of operators are considered; this is in a 
sense the case with representations whose rank is 0{N). 

More properly, we may view large N asa. thermodynamical limit [13] , since the 
number of degrees of freedom goes to infinity faster than any other physical 

parameter, and as a consequence the "macroscopic" properties of the system, 
i.e., the invariant expectation values, are fixed in spite of the great number 
of different "microscopic" realizations. This realization does not rule out the 
possibihty of searching for the so-called "master field" , that is a representa- 
tive of the equivalence class of the field configurations corresponding to the 
large-A^ vacuum, such that all invariant expectation values of the factorized 
operators can be obtained by direct substitution of the master field value into 
the definition of the operators themselves [3] . 

There has been an upsurge of interest on master fields in recent years [14,15], 
triggered by new results in non-commutative probability theory applied to the 
stochastic master field introduced in Ref. [16]. 

1.3 1/N expansion of vector models in statistical mechanics and quantum 
field theory 

The first and most successful application of the approach based on the large- 
limit and the 1/A^ expansion to field theories is the analysis of vector models 
enjoying O(A^) or SU(A^) symmetry. Actually, "vector models" is a nickname 
for a wide class of different field theories, characterized by bosonic or fermionic 
Lagrangian degrees of freedom lying in the fundamental representation of the 
symmetry group (cfr. Ref. [12] and references therein). 
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A quite general feature of these models is the possibility of expressing all 
self-interactions of the fundamental degrees of freedom by the introduction 
of a Lagrange multiplier field, a boson and a singlet of the symmetry group, 
properly coupled to the Lagrangian fields, such that the resulting effective La- 
grangian is quadratic in the A'"-component fields. One may therefore formally 
perform the Gaussian integration over these fields, obtaining a form of the 
effective action which is nonlocal, but depends only on the singlet multiplier, 
acting as a collective field; in this action N appears only as a parameter. 

The considerations developed in Subs. 1.2 make it apparent that all fluctua- 
tions of the singlet field must be suppressed in the large- limit (no residual 
degeneracy is left in the trivial representation). As a consequence, solving the 
models in this limit simply amounts to finding the singlet field configuration 
minimizing the effective action. The problem of nonlocality is easily bypassed 
by the consideration that translation invariance of the physical expectation 
values requires the action-minimizing field configuration to be invariant in 
space-time; hence the saddle-point equations of motion become coordinate- 
independent and all nonlocality disappears. 

As one may easily argue from the above considerations, the large- A" solution of 
vector models describes some kind of Gaussian field theory. Nevertheless, this 
result is not as trivial as one might imagine, since the free theory realization 
one is faced with usually enjoys quite interesting properties, in comparison 
with the naive Lagrangian free fields. Typical phenomena appearing in the 
large- A^ limit are an extension of the symmetry and spontaneous mass gen- 
eration. Moreover, when the fundamental fields possess some kind of gauge 
symmetry, one may also observe dynamical generation of propagating gauge 
degrees of freedom; this is the case with two-dimensional CP^~^ models and 
their generalizations [17,18]. 

The existence of an explicit form of the effective action offers the possibility 
of a systematic expansion in powers of 1/A^. The effective vertices of the the- 
ory turn out to be Feynman integrals over a single loop of the free massive 
propagator of the fundamental field. In two dimensions, where the physical 
properties of many vector models are especially interesting (e.g., asymptotic 
freedom), these one- loop integrals can all be computed analytically in the 
continuum version, and even on the lattice many analytical results have been 
obtained. 

The 1/A^ expansion is the starting point for a systematic computation of crit- 
ical exponents, which are nontrivial in the range 2 < d < 4, for the study of 
renormalizability of superficially nonrenormalizable theories in the same di- 
mensionality range, and for the computation of physical amplitudes. Notable 
is the case of the computation of amplitude ratios, which are independent of 
the couphng in the scaling region, and therefore are functions of 1/A" alone; 
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hopefully, their 1/N expansion possesses a nonvanishing convergence radius. 
The expansion was also useful to explore the double-scaling limit prop- 
erties of vector models [19-21]. 

The properties of the large- A?^ limit and of the 1/N expansion of continuum 
and lattice vector models were already reviewed by many authors. We there- 
fore shall not discuss this topic further. We only want to stress that this kind 
of studies can be very instructive, given the physical interest of vector mod- 
els as realistic prototypes of critical phenomena in two and three dimensions 
and as models for dynamical Higgs mechanism in four dimensions. Moreover, 
some of the dynamical properties emerging mainly from the large- studies of 
asymptotically free models (in two dimensions) may be used to mimic some of 
the features of gauge theories in four dimensions; however, at least one of the 
essential aspects of gauge theories, the presence of matrix degrees of freedom 
(fields in the adjoint representation), cannot be captured by any vector model. 

1.4 l/N expansion of matrix models: planar diagrams 

The first major result concerning the large- limit of matrix- valued field the- 
ories was due to G. 't Hooft, who made the crucial observation that, in the 
1/N expansion of continuum gauge theories, the set of Feynman diagrams con- 
tributing to any given order admits a simple topological interpretation. More 
precisely, by drawing the U(A^) fundamental fields ("quarks") as single lines 
and the U(iV) adjoint fields ("gluons") as double lines, each line carrying one 
color index, a graph corresponding to a nth-order contribution can be drawn 
on a genus n surface (i.e., a surface possessing n "holes"). In particular, the 
zeroth-order contribution, i.e., the large- A'^ limit, corresponds to the sum of all 
planar diagrams. The extension of this topological expansion to gauge mod- 
els enjoying O(A^) and Sp(2A^) symmetry has been described by Cicuta [22]. 
Large- A^ universality among O(A^), U(A^), and Sp(2A^) lattice gauge theories 
has been discussed by Lovelace [23] . 

This property has far-reaching consequences: it allows for reinterpretations 
of gauge theories as effective string theories, and it offers the possibility of 
establishing a connection between matrix models and the theory of random 
surfaces, which will be exploited in the study of the double-scaling limit. 

As a byproduct of this analysis, 't Hooft performed a summation of all planar 
diagrams in two-dimensional continuum Yang-Mills theories, and solved QCD2 
to leading nontrivial order in 1/N, finding the meson spectrum [24,25]. 

Momentum-space planarity has a coordinate-space counterpart in lattice gauge 
theories. It is actually possible to show that, within the strong-coupling ex- 
pansion approach, the planar diagrams surviving in the large- A?" limit can be 
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identified with planar surfaces built up of plaquettes by gluing them along 
half-bonds [26-28]. This construction however leads quite far away from the 
simplest model of planar random surfaces on the lattice originally proposed 
by Weingarten [29,30], and hints at some underlying structure that makes a 
trivial free-string interpretation impossible. 



1.5 The physical interpretation: QCD phenomenology 

The sum of the planar diagrams has not till now been performed in the phys- 
ically most interesting case of four-dimensional SU(A^) gauge theories. It is 
therefore strictly speaking impossible to make statements about the relevance 
of the large- limit for the description of the physically relevant case = 3. 
However, it is possible to extract from the large- analysis a number of qual- 
itative and semi-quantitative considerations leading to a very appealing pic- 
ture of the phenomenology predicted by the expansion of gauge theories. 
These predictions can be improved further by adopting Veneziano's form of 
the large- limit [31], in which not only the number of colors A^ but also the 
number of flavors Nf is set to infinity, while their ratio N/Nf is kept finite. 
We shall not enter a detailed discussion of large- A^ QCD phenomenology, but 
it is certainly useful to quote the relevant results. 



1.5.1 The large-N property of mesons 

Mesons are stable and nonintcracting; their decay amplitudes are 0{N~^^'^), 
and their scattering amplitudes are 0{N~^). 

Meson masses are finite. 

The number of mesons is infinite. 

Exotics are absent and Zweig's rule holds. 



1.5.2 The large-N property of glueballs 

Glueballs are stable and nonintcracting, and they do not mix with mesons; a 
vertex involving k glueballs and n mesons is 0{N^~^~'^/'^). 

The number of glueballs is infinite. 
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1.5.3 The large-N property of baryons 

A large- iV baryon is made out of N quarks, and therefore it possesses peculiar 
properties, similar of those of solitons [32]. 

Baryon masses are 0{N). 

The splitting of excited states is 0{1). 

Baryons interact strongly with each other; typical vertices are 0{N). 
Baryons interact with mesons with 0(1) couplings. 

1.5.4 The Tj' mass formula 

The spontaneous breaking of the SU(iVj) axial symmetry in QCD gives rise to 
the appearance of a multiplet of light pseudoscalar mesons. This symmetry- 
breaking pattern was explicitly demonstrated in the context of large- QCD 
by Coleman and Witten [33]. However, the singlet pseudoscalar is not light, 
due to the anomaly of the U(l) axial current. Since the anomaly equation 

Sat 

d,4-J^'^P^'^F"' (1-3) 

has a vanishing right-hand side in the limit Nc ^ 00 with Nf and g'^N^ fixed 
(the standard large- iV limit of non-Abelian gauge theories), the leading-order 
contribution to the mass of the rj' should be 0{1/Nc). The proportionality 
constant should be related to the symmetry-breaking term, which in turn is 
related to the so-called topological susceptibility, i.e., the vacuum expectation 
value of the square of the topological charge. The resulting relationship shows 
a rather satisfactory quantitative agreement with experimental and numerical 
results [34-38]. 

1.6 The physical interpretation: two-dimensional quantum gravity 

In the last ten years, a new interpretation of the 1/A^ expansion of matrix 
models has been put forward. Starting from the relationship between the order 
of the expansion and the topology of two-dimensional surfaces on which the 
corresponding diagrams can be drawn, several authors [39-43] proposed that 
large- iV matrix models could provide a representation of random lattice two- 
dimensional surfaces, and in turn this should correspond to a realization of 
two-dimensional quantum gravity. These results were found consistent with 
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independent approaches, and proper modifications of the matrix self-couphngs 
could account for the incorporation of matter. 

The functional integrals over two-dimensional closed Riemann manifolds can 
be replaced by the discrete sum over all (piecewise flat) manifolds associated 
with triangulations. It is then possible to identify the resulting partition func- 
tion with the vacuum energy 

E^^-logZ^, (1.4) 



obtained from a properly defined N xN matrix model, and the topological ex- 
pansion of two-dimensional quantum gravity is nothing but the 1/N expansion 
of the matrix model. 

The partition function of two-dimensional quantum gravity is expected to 
possess well-defined scaling properties [44]. These may be recovered in the 
matrix model by performing the so-called "double-scahng hmit" [45-47] . This 
limit is characterized by the simultaneous conditions 

TV^oo, g-^Qc, (1-5) 



where g is a, typical self-coupling and g^ is the location of some large- phase 
transition. The limits are however not independent. In order to get nontrivial 
results, one is bound to tune the two conditions (1.5) in such a way that the 
combination 

X ^ {g - g,)N^''^' (1.6) 



is kept finite and fixed. 71 is a computable critical exponent, usually called 
"string susceptibility" . According to Ref. [44], it is related to the central charge 
c of the model by 
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25 - c + ^(l-c)(25-c) 



(1.7) 



An interesting rcintcrprctation of the double-scaling limit relates it to some 
kind of finite-size scaling in a space where N plays the role of the physical 
dimension L [21,48,49]. Research in this field has exploded in many directions. 
A wide review reflecting the state of the art as of the year 1993 appeared in 
the already- mentioned volume by Brezin and Wadia [8]. Here wc shall only 
consider those results that are relevant to our more restricted subject. 
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2 Unitary matrices 



2. 1 General features of unitary-matrix models 

Under the header of unitary-matrix models we class all the systems character- 
ized by dynamical degrees of freedom that may be expressed in terms of the 
matrix representations of the unitary groups U(A'^) or special unitary groups 
SU(A^) and by interactions enjoying a global or local U(A^)l x V{N)ii sym- 
metry. Typically we shall consider lattice models, with no restriction on the 
lattice structure and on the number of lattice points, ranging from 1 (single- 
matrix problems) to infinity (infinite-volume limit) in an arbitrary number of 
dimensions. 

In the field-theoretical interpretation, i.e., when considering models in infinite 
volume and in proximity of a fixed point of some (properly defined) renormal- 
ization group transformation, such models will have a continuum counterpart, 
which in turn shall involve unitary-matrix valued fields in the case of spin 
models, while for gauge models the natural continuum representation will be 
in terms of hermitian matrix (gauge) fields. 

A common feature of all unitary-matrix models will be the group-theoretical 
properties of the functional integration measure: for each dynamical variable 
the natural integration procedure is based on the left- and right-invariant Haar 
measure 

d//(C/) = dii{UV) = d/i(yC/), j dii(U) = 1. (2.1) 

An explicit use of the invariance properties of the measure and of the interac- 
tions (gauge fixing) can sometimes lead to formulations of the models where 
some of the symmetries are not apparent. Global U(A^) invariance is how- 
ever always assumed, and the interactions, as well as all physically interesting 
observables, may be expressed in terms of invariant functions. 

It is convenient to introduce some definitions and notations. An arbitrary 
matrix representation of the unitary group U(A^) is denoted by T)^^^(U). The 
characters and dimensions of irreducible representations are X{r){U) = T>^^}{U) 
and d{r) respectively, (r) is characterized by two set of decreasing positive 
integers {/} = h, ...Is and {m} = mi, rrit. We may define the ordered set of 
integers {A} = Ai,...,Aivby the relationships 

Afc = /fe, {k=l,...,s), Afe = 0, {k = s + l,...,N -t), 

Afe = -mjv-fc+i, {k^N -t + l,...,N). (2.2) 
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It is then possible to write down explicit expressions for all characters and 
dimensions, once the eigenvalues exp i^j of the matrix U are known: 



det||exp{i0^(Aj + A^-i)}|| 
= det||exp{i0.(iV-,-)}|| ' ^2.3) 

^(.)- -X(A)(1)- (2.4) 



The general form of the orthogonality relations is 

/ MU)V^i{U)V^:!*{U) = -]-5r,sKcKd. (2.5) 



Further relations can be found in Ref . [50] . 

The matrix Uab itself coincides with the fundamental representation (1) of the 
group, and enjoys the properties 

X(i){U) = Tr [/, = TV, UabUl = 5,, . (2.6) 



The measure d/i([/) (which we shall also denote simply by dU), when the 
integrand depends only on invariant combinations, may be expressed in terms 
of the eigenvalues [51]. 



2.2 Chiral models and lattice gauge theories 



Unitary matrix models defined on a lattice can be divided into two major 
groups, according to the geometric and algebraic properties of the dynamical 
variables: when the fields are defined in association with lattice sites, and the 
symmetry group is global, i.e., a single U{N)l x U{N)r transformation is 
applied to all fields, we are considering a spin model (principal chiral model); 
in turn, when the dynamical variables are defined on the links of the lattice 
and the symmetry is local, i.e., a different transformation for each site of the 
lattice may be performed, we are dealing with a gauge model (lattice gauge 
theory). As we shall see, these two classes are not unrelated to each other: 
an analogy between d-dimensional chiral models and 2(i-dimensional gauge 
theories can be found according to the following correspondence table [52] : 
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spin 



gauge 



site, link 



link, plaquette 



loop 



surface 



length 



area 



mass 



string tension 



two-point correlation Wilson loop 

While this correspondence in arbitrary dimensions is by no means rigorous, 
there is some evidence supporting the analogy. 

In the case d — 1, which we shall carefully discuss later, one can prove an 
identity between the partition function (and appropriate correlation functions) 
of the two-dimensional lattice gauge theory and the corresponding quantities of 
the one-dimensional principal chiral model. Both theories are exactly solvable, 
both on the lattice and in the continuum limit, and the correspondence can 
be explicitly shown. 

Approximate real-space renormalization recursion relations obtained by Migdal 
[53] are identical for ci-dimensional chiral models and 2(i-dimensional gauge 
models. 

The two-dimensional chiral model and the (phenomenologically interesting) 
four- dimensional non-Abelian gauge theory share the property of asymptotic 
freedom and dynamical generation of a mass scale. In both models these prop- 
erties are absent in the Abehan case (XY model and U(l) gauge theory re- 
spectively), which shows no coupling-constant renormalization in perturbation 
theory. 

The structure of the high-temperature expansion and of the Schwinger-Dyson 
equations is quite similar in the two models. 

It will be especially interesting for our purposes to investigate the Schwinger- 
Dyson equations of unitary-matrix models and discuss the peculiar properties 
of their large- limit. 

2.3 Schwinger-Dyson equations in the large-N limit 

In order to make our analysis more concrete, we must at this stage consider 
specific forms of interactions among unitary matrices, both in the spin and in 
the gauge models. The most dramatic restriction that we are going to impose 
on the lattice action is the condition of considering only nearest-neighbor 
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interactions. The origin of this restriction is mainly practical, because non 
nearest-neighbor interactions lead to less tractable problems. We assume that, 
for the systems we are interested in, it will always be possible to find a lattice 
representation in terms of nearest-neighbor interactions within the universality 
class. 

Let us denote by x an arbitrary lattice site, and by x, fi an arbitrary lattice link 
originating in the site x and ending in the site x + /i: /i is one of the d positive 
directions in a d-dimensional hypercubic lattice. A plaquette is identified by 
the label x, /i, v, where the directions /x and v {[x ^ v) specify the plane 
where the plaquette lies. The dynamical variables (which we label by U in the 
general case) are site variables in spin models and link variables U^^^i in 
gauge models. 

The general expression for the partition function is 

Z = jjlMU)^M-(^S{U)i (2.7) 

where j3 is the inverse temperature (inverse coupling) and the integration is 
extended to all dynamical variables. The action S{U) must be a function en- 
joying the property of extensivity and of (global and local) group invariance, 
and respect the symmetry of the lattice. Adding the requisite that the inter- 
actions involve only nearest neighbors, we find that a generic contribution to 
the action of spin models must be proportional to 

Exw(t/.^i+^)+h.c., (2.8) 

and for gauge models to 

E X{r){U,,^U,+^,,uU,,,Ul,) + h.c. , (2.9) 

X,IJ.,V 

where (r) is in principle arbitrary, and the summation is extended to all ori- 
ented links of the lattice in the spin case, to all the oriented plaquettes in the 
gauge case. In practice we shall mostly focus on the simplest possible choice, 
corresponding to the fundamental representation. In order to reflect the ex- 
tensivity of the action, i.e., the proportionality to the number of space and 
internal degrees of freedom, it will be convenient to adopt the normalizations 

S{U) = - E ^(Tr U^UU^ + h.c.) (spin), (2.10) 
S{U) = - E ^(Tr U,,,U,+,,Mi+.,^Ul, + h.c.) (gauge). (2.11) 

x,ii,u 
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Once the lattice action is fixed, it is easy to obtain sets of Schwinger-Dyson 
equations relating the correlation functions of the models. These are the quan- 
tum field equations and solving them corresponds to finding a complete solu- 
tion of a model. It is extremely important to notice the simplifications occur- 
ring in the Schwinger-Dyson equations when the large- A?" limit is considered. 
These simplifications are such to allow, in selected cases, explicit solutions to 
the equations. 

Before proceeding to a derivation of the equations, we must preliminarily 
identify the sets of correlation functions we are interested in. For obvious 

reasons, these correlations must involve the dynamical fields at arbitrary space 
distances, and must be invariant under the symmetry group of the model. 
Without pretending to achieve full generality, we may restrict our attention 
to such typical objects as the invariant correlation functions of a spin model 

= ^ (Trf[U.^Ul}i (2.12) 

and to the so-called Wilson loops of a gauge model 

iy(C) = i^^TVn^z), (2.13) 

where C is a closed arbitrary walk on the lattice, and Yliec ^^e ordered 
product over all the links along the walk. It is worth stressing that the action 
itself is a sum of elementary Green's functions (elementary Wilson loops). 

More general invariant correlation functions may involve expectation values of 
products of invariant operators similar to those appearing in the r.h.s. of Eqs. 
(2.12) and (2.13). The already mentioned property of factorization allows us 
to express the large- A?" limit expectation value of such products as a product of 
expectation values of the individual operators. As a consequence, the large- A?" 
form of the Schwinger-Dyson equations is a (generally infinite) set of equations 
involving only the above-defined quantities. 

For sake of clarity and completeness, we present the exphcit large- A^ form 
of the Schwinger-Dyson equations for the models described by the standard 
actions (2.10) and (2.11). For principal chiral models [54], 

= G(")(xi,yi, ...,Xn,yn) 
-F/?^[G^"+^^(a;i,Xi ...,x„,y„) - G^'"\xi + ii,yi, ...,x„,y„) 

n 

+ Yl[^xuxs G^'~^\xi, yi, Xs-l, ys-i) G^''~'+^\xs, ys, ■-, Xn, yn) 

s=2 
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- ^xi,ys G^^\xi,yi, Xs, Vs) C^" ^\xs+l, Vs+l: ■-, ^^n, Vn)] ■ (2-14) 

For lattice gauge theories [55,56], 



P[i:WiC.,,u)-WiC,^,,,^)] = Y.S,,yWiC,,y)WiCy,,), (2.15) 

where W{Cx,n„) is obtained by replacing C4,^ with Ux,^iUa:+^^Ml+,y,^ in the loop 
C, and Cy^x are the sub-loops obtained by splitting C at the intersection 
point, including the "trivial" splitting. Eqs. (2.15) are commonly known as the 
lattice Migdal-Makeenko equations. The derivation of the Schwinger-Dyson 
equations is obtained by performing infinitesimal variations of the integrand 
in the functional integral representation of expectation values and exploiting 
invariance of the measure. 



2.4 Survey of different approaches 



Schwinger-Dyson equations are the starting point for most techniques aiming 
at the explicit evaluation of large- A'^ vacuum expectation values for nontrivial 
unitary-matrix models. The form exhibited in Eqs. (2.14) and (2.15) involves 
in principle an infinite set of variables, and it is therefore not immediately 
useful to the purpose of finding explicit solutions. 

Successful attempts to solve large- A?" matrix systems have in general been 
based on finding reformulations of Schwinger-Dyson equations involving more 

restricted sets of variables and more compact representations (collective fields) . 
As a matter of fact, in most cases it turned out to be convenient to define gen- 
erating functions, whose moments are the correlations we are interested in, and 
whose properties are usually related to those of the eigenvalue distributions 
for properly chosen covariant combinations of matrix fields. 

By "covariant combination" we mean a matrix-valued variable whose eigen- 
values are left invariant under a general SU(A^) x SU(A^) transformation of 
the Lagrangian fields. Such objects are typically those appearing in the r.h.s. 
of Eqs. (2.12) and (2.13) before the trace operation is performed. Under the 
SU(A^) X SU(A^) transformation U — > VUW'^ . these operators transform ac- 
cordingly to C — > VOV\ and therefore their eigenvalue spectrum is left un- 
changed. 

Without belaboring on the details (some of which will however be exhibited in 
the discussion of the single-link integral presented in Sect. 3), we only want to 
mention that the approach based on extracting appropriate Schwinger-Dyson 
equations for the generating functions is essentially algebraic in nature, involv- 



17 



ing weighted sums of infinite sets of equations in the form (2.14) or (2.15), 
identification of the relevant functions, and resohition of the resulting alge- 
braic equations, where usually a number of free parameters appear, whose 
values are fixed by boundary and/or asymptotic conditions and analyticity 
constraints. The approach based on direct replacement of the eigenvalue dis- 
tributions in the functional integral and the minimization of the resulting 
effective action leads in turn to integral equations which may be solved by 
more or less straightforward techniques. These two approaches are however 
intimately related, since the eigenvalue density is usually connected with the 
discontinuity along some cut in the complex-plane extension of the generating 
function, and one may easily establish a step-by-step correspondence between 
the algebraic and functional approach. 

Let us finally mention that the procedure based on introducing invariant de- 
grees of freedom and eigenvalue density operators has been formalized by 
Jevicki and Sakita [57,58] in terms of a "quantum collective field theory", 
whose equations of motion are the Schwinger- Dyson equations relevant to the 
problem at hand. 

A quite different application of the Schwinger-Dyson equations is based on the 
strong-coupling properties of the correlation functions. In the strong-coupling 
domain, expectation values are usually analytic in the coupling f3 within some 
positive convergence radius, and their boundary value at = can easily be 
evaluated. As a consequence, it is formally possible to solve Eqs. (2.14) and 
(2.15) in terms of strong-coupling series by sheer iteration of the equations. 
This procedure may in practice turn out to be too cumbersome for practical 
purposes; however, in some circumstances, it may lead to rather good approxi- 
mations [59,60] and even to a complete strong-coupling solution. Continuation 
to the weak-coupling domain is however a rather nontrivial task. 

As a special application of the strong-coupling approach, we must mention the 
attempt (pioneered by Kazakov, Kozhamkulov and Migdal [61]) to construct 
an effective action for the invariant degrees of freedom by means of a modified 
strong-coupling expansion, and explore the weak-coupling regime by solving 
the saddle-point equations of the resulting action. This technique might be 
successful at least in predicting the location and features of the large- phase 
transition which is relevant to many physical problems, as mentioned in Sect. 
1. 

A numerical approach to large- A" lattice Schwinger-Dyson equations based on 
the minimization of an effective large- A^ Fokker-Plank potential and suited for 
the weak-coupling regime was proposed by Rodrigues [62] . 

Another relevant application of the Schwinger-Dyson equations is found in the 
realm of the so-called "reduced" models. These models, whose prototype is the 
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Eguchi-Kawai formulation of strong-coupling large- lattice gauge theories 
[63], are based on the physical intuition that, in the absence of fluctuations, due 
to translation invariance, the space extension of the lattice must be essentially 
irrelevant in the large- A'^ limit, since all invariant physics must be already 
contained in the expectation values of (properly chosen) purely local variables. 
More precisely, one might say that, when iV — > cxo, the SU(A'^) group becomes 
so large that it accommodates the full Poincarc group as a subgroup, and in 
particular it should be possible to find representations of the translation and 
rotation operators among the elements of SU(A^). As a consequence, one must 
be able to reformulate the full theory in terms of a finite number of matrix 
field variables defined at a single space-time site (or on the d links emerging 
from the site in the case of a lattice gauge theory) and of the above-mentioned 
representations of the translation group. This reformulation is called "twisted 
Eguchi-Kawai" reduced version of the theory [64,65]. 

We shall spend a few more words on the reduced models in Sect. 7. Moreover, a 
very good review of their properties has already appeared many years ago [6]. 
In this context, we must only mention that the actual check of validity of the 
reduction procedure is based on deriving the Schwingcr-Dyson equations of the 
reduced model and comparing them with the Schwinger-Dyson equations of 
the original model. Usually the equivalence is apparent already at a superficial 
level when naively applying to correlation functions of the reduced model 
the symmetry properties of the action itself. This procedure however requires 
some attention, since the limit of infinitely many degrees of freedom within 
the group itself allows the possibility of spontaneous breakdown of some of 
the symmetries which would be preserved for any finite value of N. In this 
context, we recall once more that large A^ is a thermodynamical limit: A^ must 
go to infinity before any other limit is considered, and sometimes the limiting 
procedures do not commute. It is trivial to recognize that, when the strong- 
coupling phase is considered, symmetries are unbroken, and the equivalence 
between original and reduced model may be established without further ado. 
Problems may occur in the weak-coupling side of a large- A?" phase transition. 

An unrelated and essentially numeric approach to solving the large- A^ limit of 
lattice matrix models is the coherent state variational algorithm introduced by 
Yaffe and coworkers [66,67]. We refer to the original papers for a presentation 
of the results that may be obtained by this approach. 



3 The single-link integral 
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3.1 The single-link integral in external field: finite-N solution 



All exact and approximate methods of evaluation of the functional integrals 
related to unitary-matrix models must in principle face the problem of per- 
forming the simplest of all relevant integrations: the single-link integral. The 
utmost importance of such an evaluation makes it proper to devote to it an 
extended discussion, which will also give us the opportunity of discussing in a 
prototype example the different techniques that may be applied to the models 
we are interested in. 

A quite general class of single-link integrals may be introduced by defining 

= y dC/ ex.p[N Tr(AtC/ + C/t A)] , (3.1) 

where as usual U is an element of the group U(A^) and A is now an arbitrary 
N X N matrix. The U(A'") invariance of the Haar measure implies that the 
one link integral (3.1) must depend only on the eigenvalues of the Hermitian 
matrix A'' A, which we shall denote by Xi,...,Xn- The function Z{xi, ...,Xn) 
must satisfy a Schwinger-Dyson equation: restricting the variables to the U(A^) 
singlet subspace, the Schwinger-Dyson equation was shown to be equivalent 
to the partial differential equation [68,69] 



J_ ^dZ_ 1 ^ Xs (dZ dZ\_^ 

N'^ dxl ^ 9xk N"^ f^k^k- Xs \dxk dxs) 



(3.2) 



with the boundary condition Z(0, ...,0) = 1 and the request that Z be com- 
pletely symmetric under exchange of the Xj. 

It is convenient to reformulate the equation in terms of the new variables 
Zk — 2Ny/xk, and to parameterize the solution in terms of the completely 
antisymmetric function Z{zi, ...,zn) by defining 



The equation satisfied by Z can be shown to reduce to 



■ k k k k k 



Z^O. 
(3.4) 
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Eq. (3.4) has the structure of a fermionic many-body Schrodinger equation. 
With some ingenuity it may be solved in the form of a Slater determinant of 
fermion wavefunctions. In conclusion, we obtain, after proper renormalization 
[70] (see also [71]), 

Z(.„ z^) = 2^(^-)/^ fn k^) "^^"^^t'^fK (3.5) 

\fe=o / detll^;^-^ ^11 

where Ii{z) is the modified Bessel function. Eq. (3.5) is therefore a representa- 
tion of the single-link integral in external field for arbitrary U(A^) groups. By 
taking proper derivatives with respect to its arguments one may in principle 
reconstruct all the cumulants for the group integration of an arbitrary string 
of (uncontracted) matrices [72,73]. 

Some special limits of the general expression (3.5) may prove useful. Let us 
first of all consider the case when A is proportional to the identity matrix: 
A = al and therefore Zi = 2Na and 

Z(2Na, ...,2Na) = det \\Ii^j{2Na)\\. (3.6) 

As we shall see, this is exactly Bars' and Green's solution for U(A^) lattice 
gauge theory in two dimensions [74] . 

When only one eigenvalue of A is different from zero the result is 

Z{2Na, 0, 0) = (A^ - 1)! (Na)^'^ lN-i(2Na). (3.7) 

The large- limit will be discussed in the next subsection. 
3.2 The external field problem: large-N limit 

For our purposes it is extremely important to extract the limiting form of 
Eq. (3.5) when A^ — > oo. In principle, it is a very involved problem, since the 
dependence on A^ comes not only through the Zi but also from the dimension 
of the matrices whose determinant we must evaluate. It is however possible to 
obtain the limit, either by solving separately the large- A^ version of Eq. (3.2), 
or by directly manipulating Eq. (3.5). 

In the first approach, we introduce the large- A" parameterization 

Z^expNW, (3.8) 
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where W is now proportional to A^; we then obtain from Eq. (3.2), dropping 
second-derivative terms that are manifestly depressed in the large- A?" limit [69], 



'dWV dW 1 ^ Xs fdW _ dW 
^ dxk N j^j^Xg- Xk V dxs dxk ^ 



dxh 



= 1. 



(3.9) 



It is possible to show that in the large-A^ limit Eq. (3.9) admits solutions, 
which can be parameterized by the expression 



dW 



dxk \/xk + c 



J-V I 

2N V \/^feTc + V^TTc 



c > 0. (3.10) 



Substitution of Eq. (3.10) into Eq. (3.9) and some algebraic manipulation lead 
to the consistency condition 



2Ar^ v^TT^ 



0, 



(3.11) 



which in turn admits two possible solutions: 
a) c is determined by the condition 



(3.12) 



implying c < |; this is a "strong coupling" phase, requiring that the eigenval- 
ues satisfy the bound 



i.e., at least some of the Xs are sufficiently small; 
h) when 



2N^ 



(3.14) 



then the solution corresponds to the choice c = 0; this is a "weak coupling" 
phase, and all eigenvalues are large enough. 

Direct integration of Eq. (3.10) with proper boundary conditions leads to the 
large- result [69] 
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= 2 ^ - — ^ log(x/i^ + V^7+^) - ATc - - A^, 

(3.15) 

which must be supplemented with Eq. (3.12) in the strong- coupling regime 
(3.13), while c — reproduces the weak-coupling result by Brower and Nauen- 
berg. Amazingly enough, setting c = in Eq. (3.15) one obtains the naive 
one-loop estimate of the functional integral, which turns out to be exact in 
this specific instance. 

It is possible to check that Eq. (3.15) is reproduced by carefully taking the 
large- iV limit of Eq. (3.5), which requires use of the following asymptotic limits 
of Bessel functions [70] 



l-k 



-1/4 



1 + 



h{z) 



X exp (^V + — (strong coupling) , 

1 



V27ri 



exp z (weak coupling) . 



(3.16) 
(3.17) 



An essential feature of Eq. (3.15) is the appearance of two different phases 
in the large-A^ limit of the single-link integral. Such a transition would be 
mathematically impossible for any finite value of N; however it affects the 
large- A?" behavior of all unitary-matrix models and gives rise to a number of 
interesting phenomena. A straightforward analysis of Eq. (3.15) shows that 
the transition point corresponds to the condition 

*^wSi= = l. (3.18) 



It is also possible to evaluate the difference between the strong- and weak- 
coupling phases of W in the neighborhood of t = 1, finding the relationship 
[69] 

W^strong - W^weaic ~ (^ - 1)'- (3.19) 



As a consequence, we may classify this phenomenon as a "third order phase 
transition" . 



23 



3.3 The properties of the determinant 



The large- factorization of invariant amplitudes is a well-established prop- 
erty of products of operators defined starting from the fundamental represen- 
tation of the symmetry group. Operators corresponding to highly nontrivial 
representations may show a more involved pattern of behavior in the large- 
limit. Especially relevant from this point of view are the properties of de- 
terminants of covariant combinations of fields [52,75]; we will consider the 
quantities 

A(a;) = det \UoUl\ (3.20) 

for lattice chiral models and 

A(C) = det n Ui (3.21) 

i&c 

for lattice gauge theories. 

The expectation values of these operators may act as an order parameter for 
the large- phase transition characterizing the class of models we arc taking 
into consideration. Indeed the determinant picks up the phase characterizing 
the U(l) subgroup that constitutes the center of U(A^). Moreover, since 

U(iV).U(l)x^, 

SU(iV) — > U(A^) as iV — > cxo because Z^v — >■ U(l); therefore the determinant 
of the U(A^) theory in the large- hmit reflects properties of the center of 
SU(Ar). 

In lattice models this Abelian U(l) subgroup is not decoupled, as it happens 
in the continuum theory, and therefore (A) does not in general have on the 
lattice the free-theory behavior it has in the continuum. 

The basic properties of the determinant may be explored by focusing once 
more on the external field problem we discussed above. Let us introduce a 
class of determinant operators, and define their expectation values as [76] 

A(0 - /det - /d^dct^'exp[ArTr(mA + At^)] 

In order to parameterize the SU(A^) external-source integral, besides the eigen- 
values Xi of AA^ , a new external parameter must be introduced, that couples 



24 



to the determinant: 



^ = -i- (log det - log det A) . 



(3.23) 



Because of the symmetry properties, A^'^ may only depend on the eigenvalues 
z and on 9. It was found that, when U enjoys U(A^) symmetry (with finite N), 



=exp(iM^)^, 



(3.24) 



where Zi is the solution of the following Schwinger- Dyson equation, general- 
izing Eq. (3.4): 



1 

TV 



E^^^ + (3-2Ar)^,,^-^,^ 

k k /c 



+ -N{N-l){N-2) 



(3.25) 



Zi satisfy the property 



zi = mzA n 



1 d 



Zq = det \\zl ^/j_i_/(z 



(3.26) 



When the weak-coupling condition t = J2k ^/^k < 1 is satisfied, the leading 
contribution to the large- A?^ limit of all Zi is the same: 



z, 



exp 



Y^^k- 27r2;fc + log{zi - Zk) 

k k i<k 



(3.27) 



In order to determine the large- A?" limit of A^'^, one therefore needs to com- 
pute the 0(1) factor in front of the exponentially growing term (3.27). It is 
convenient to define 



(3.28) 



whose Schwinger-Dyson equation may be extracted from Eq. (3.25) and takes 
the form 



25 



1 

N 
/2 



k 



ZkZi 



dXi dXi 



dzi 



(3.29) 



Let us introduce the large- A?" Ansatz 



(3.30) 



reducing Eq. (3.29) to 



1 ^ 1 A^Xi , ,dXz 



I' -AX, 



(3.31) 



Removing terms that are depressed by two powers of l/A'^, we are left with a 
consistent equation whose solution is 



X, = (l-t)i('M). 



(3.32) 



Finally we can compute the weak-coupling large- limit of A*^'); 

A(') — > exp{iNW) (1 -t)^^\ t< 1. 

AT— ►oo 



(3.33) 



Prom the standard strong-coupling expansion we may show that 



A(^) — > when t > 1. 

AT— >oo 



(3.34) 



An explicit evaluation, starting from the exact expression (3.26), expanded in 

powers of 1/zk for arbitrary N, allows us to show that the quantities Zi may 
be obtained from Eqs. (3.27) and (3.28) by expanding Eq. (3.32) up to 2nd 
order in t with no 0(l/iV^) corrections. A*^''' according to this result violate 
factorization; in turn, they take the value which would be predicted by an 
effective Gaussian theory governing the U(l) phase of the field U. 



3.4 Applications to mean field and strong coupling 



The single-link external-field integral has a natural domain of application in 
two important methods of investigation of lattice field theories: mean-field 
and strong-coupling expansion. Extended papers and review articles have been 
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devoted in the past to these topics (cfr. Ref. [77] and references therein), and 
we shall therefore focus only on those results that are specific to the large- A?^ 
limit and to the expansion. 

Let us first address the issue of the mean-field analysis, considering for sake 
of definiteness the case of rf-dimensional chiral models, but keeping in mind 
that most results can be generalized in an essentially straightforward manner 
to lattice gauge theories. The starting point of the mean-field technique is the 
application of the random field transform to the functional integral: 



Zjv = / dUn expj AT/? g Tr {UnUl+^ + U^+M) } 
= / dK.d^„exp|iV/?^Tr(KK!+M + ^n+^KJ) 
-7V^Tr(AKJ + K^9| 

n J 

X j dU^ exp jiV ^ Tr (A^UI + t/^^t ) | , (3.35) 

where Vn and are arbitrary complex NxN matrices. Therefore the in- 
tegration over Un is just the single-link integral we discussed above. As a 
consequence, the original chiral model is formally equivalent to a theory of 
complex matrices with effective action 



- ^ s^^A, y) = /3 E ^ (kkVm + k+mK!) - E 1^ + k4) 

+ Y.mAnAi). (3.36) 

n 

The leading order in the mean-field approximation is obtained by applying 
saddle-point techniques to the effective action, assuming saddle-point values 
of the fields A„ and Vn that are translation-invariant and proportional to the 
identity. 

We mention that, in the case at hand, the large- saddle-point equations in 
the weak-coupling phase are: 

A„ = a = 2pdv, Vn^v^l - ^, (3.37) 
and they are solved by the saddle-point values 
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leading to a value of the free and internal energy 

Fd 1, ^ 1 1 d Fd 



a--log2a--, (3.39) 



N^L 2 2' 2ddpN^L 



The strong-coupling solution is trivial: v — a — 0, and there is a first-order 
transition point at 

Pcd^l Vc^i CLc^i (3.40) 



One may also compute the quadratic fluctuations around the mean-field saddle 
point by performing a Gaussian integral, whose quadratic form is related to 
the matrix of the second derivatives of W with respect to the fields, and 
generate a systematic loop expansion in the effective action (3.36), which in 
turns appears to be ordered in powers oil/d. Therefore mean-field methods are 
especially appropriate for the discussion of models in large space dimensions, 
and not very powerful in the analysis oi d — 2 models. The very nature of 
the transition cannot be taken for granted, especially at large A^. However, 
when d > 3 there is independent evidence of a first-order phase transition for 
> 3. We mention that a detailed mean- field study of SU(iV) chiral models 
in d dimensions appeared in Refs. [78,79]. 

When willing to extend the mean-field approach, it is in general necessary to 
find a systematic expansion of the functional 1^(74^4^) in the powers of the 
fluctuations around the saddle-point configurations. Moreover, one may choose 
to consider not only the large- value of the functional, but also its expansion 
in powers if in order to make predictions for large but flnite values of 

N. The expansion of Wq up to fourth order in the fluctuations was performed 
in Ref. [80], where explicit analytic results can be found. A technique for the 
weak-coupling 1/A^^ expansion of W can be found in Ref. [81]. We quote the 
complete 0{1/N'^) result: 



W_J_ 



E - ^ E log ^ --N^ + iog(i - t)-v« 

. 1 a,b 



a 



(3.41) 



where t = J2a V-^a- ^'i- (3-41) can also be expanded in the fluctuations around 
a saddle-point configuration. Extension to SU(A^) with large A" was also con- 
sidered. A discussion of large- A^ mean field for lattice gauge theories can be 
found in Refs. [79,82-85]. 

Let us now turn to a discussion of the main features of the large- A^ strong- 
coupling expansion. A preliminary consideration concerns the fact that it is 
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most convenient to reformulate the strong-coupling expansion (i.e., the expan- 
sion in powers of /5) into a character expansion, which is ordered in the number 
of lattice steps involved in the effective path that can be associated with each 
nontrivial contribution to the functional integral. The large- character ex- 
pansion will be discussed in greater detail in Subs. 4.8. Here we only want to 
discuss those features that are common to any attempt aimed at evaluating 
strong-coupling series for expectation values of invariant operators in the con- 
text of U(A^) and SU(A^) matrix models, with special focus on the large-iV 
behavior of such series. 

The basic ingredient of strong-couphng computations is the knowledge of the 
cumulants, i.e., the connected contributions obtained performing the invariant 
group integration of a string of uncontracted U and U"^ matrices. U(A^) group 
invariance insures us that these group integrals can be non-zero only if the 
same number of U and U"^ matrices appear in the integrand. SU(A^) is slightly 
different in this respect, and its peculiarities will be discussed later and are 
not relevant to the present analysis. 

It was observed a long time ago that the cumulants, whose group structure 
is that of invariant tensors with the proper number of indices, involve A^- 
dependent numerical coefficients. The asymptotic behavior of these coefficients 
in the large- limit was studied first by Weingarten [86]. However, for finite N, 
the coefficients written as function of N are formally plagued by the so-called 
DeWit-'t Hooft poles [87], that are singularities occurring for integer values of 
N. The highest singular value of N grows with the number n of U matrices 
involved in the integration, and therefore for sufficiently high orders of the 
series it will reach any given finite value. A complete description of the pole 
structure was presented in Ref. [72]; not only single poles, but also arbitrary 
high-order poles appear for large enough n, and analyticity is restricted to 
N > n. Obviously, since group integrals are well defined for all n and A^, this 
is only a pathology of the 1/A^ expansion. Finite- A^ results are finite, but they 
cannot be obtained as a continuation of a large- A^ strong-coupling expansion. 
However, it is possible to show that the strict N ^ oo limit of the series exists, 
and moreover, for sufficiently small 3 and sufficiently large A^, the limiting 
series is a reasonable approximation to the true result, all nonanalytic effects 
being 0{(3'^^) in U(A^) models and in SU(A^) models. As a consequence, 

computing the large- A?^ limit of the strong-coupling series is meaningful and 
useful in order to achieve a picture of the large- A^ strong-coupling behavior of 
matrix models, but the evaluation of 0{1/N'^) or higher-order corrections in 
the strong-coupling phase is essentially pointless. 

The large- A^ limit of the external-field single-link integral has been considered 
in detail from the point of view of the strong-coupling expansion. In particular, 
one may obtain expressions for the coefficients of the expansion of W in powers 
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of the moments of AA^: setting 

1 oo 

p„ = - IV(AAt)", W = Y. E (3.42) 

n=l ai,...,a„ 
J2k ^"fc=" 



one gets 



W, 



ai, 



..n C^n + Ek^^k-Sy. ^l (2fc)! 
^ (2n)! \n (^!)^ 



1 



(3.43) 



Further properties of this expansion can be found in the original reference [88] . 

A character-expansion representation of the single-hnk integral was also pro- 
duced for arbitrary U(A^) integrals in Ref. [73]. Strong-coupling expansions for 
large- A?" lattice gauge theories have been analyzed in detail by Kazakov [26,89], 
O'Brien and Zubcr [27], and Rostov [28], who proposed reinterpretations in 
terms of special string theories. 



3.5 The single-link integral in the adjoint representation 



The integral introduced at the beginning of Sect. 3 is by no means the most 
general single-link integral one can meet in unitary-matrix models. As men- 
tioned in Sect. 2, any invariant function of the C/'s is in principle a candidate 
for a lattice action. In practice, the only case that has been considered till now 
that cannot be reduced to Eq. (3.1) is the integral introduced by Itzykson and 
Zuber [50] 

/(Mi,M2)= f dUexpTT{MiUM2U^), (3.44) 



where Mi and M2 arc arbitrary Hermitian matrices. This is a special instance 
of the single-link integral for the coupling of the adjoint representation of U 
to an external field. 

The result, because of U(A^) invariance, can only depend on the eigenvalues 
mii and of the Hermitian matrices. Several authors [50,90,91] have inde- 
pendently shown that 

/(M.,M.) = (n'p!) ., (3.46) 

Vp=i / ^{'^u, ■-, miN) A(m2i, m2N) 
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where A{mi, ...,171^) = nj>j("^j ~ ""^j) is the Vandemonde determinant. A 
series expansion for I {Mi, M2) in terms of the characters of the unitary group 
takes the form 

/(Mi,M2) = ^^^XM(Mi)x(r)(M2), (3.46) 

where (7(r) is the dimension of the representation (r) of the permutation group; 
wc will present an explicit evaluation of a(^r) in Eq- (4.100). Eq. (3.45) plays a 
fundamental role in the decoupling of the "angular" degrees of freedom when 
models involving complex Hermitian matrices are considered. 

An interesting development based on the use of Eq. (3.45) is the so-called 
"induced QCD" program, aimed at recovering continuum large-A^ QCD by 
taking proper limits in the parameter space of the lattice Kazakov-Migdal 
model [92] 

S = AT^IY V^($,) - Ar^lV($,C/,,^$,+^C/tj, (3.47) 

X X,IJ, 

where Ux,ij, is the non-Abelian gauge field and is a Hermitian N x N 
(matrix- valued) Lorentz-scalar field. The Itzykson-Zuber integration (3.44) 
allows the elimination of the gauge degrees of freedom and reduces the problem 
to studying the interactions of Hermitian matrix fields (with self-interactions 
governed by the potential V). Discussion of the various related developments 
is beyond the scope of the present report. It will be enough to say that, while 
one may come to the conclusion that this model does not induce QCD, it 
is certainly related to some very interesting (and sometimes solvable) matrix 
models (cfr. Ref. [93] for a review). 



4 Two-dimensional lattice Yang-Mills theory 

4-1 Two-dimensional Yang-Mills theory as a single-link integral 

The results presented in the previous section allow us to analyze the simplest 
physical system described by a unitary-matrix model. As we shall see, one of 
the avatars of this system is a Yang- Mills theory in two dimensions (YM2), in 
the lattice Wilson formulation. Notwithstanding the enormous simplifications 
occurring in this model with respect to full QCD, still some nontrivial features 
are retained, and even in the large- limit some interesting physical properties 
emerge. It is therefore worth presenting a detailed discussion of this system. 
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which also offers the possibihty of comparing the different technical approaches 
to the large- solution in a completely controlled situation. 

The lattice formulation of the two-dimensional U(A^) gauge theory is based 
on dynamical variables [4,^* which are defined on links; however, because of 
gauge invariancc, in two dimensions there are no transverse gauge degrees of 
freedom, and a onc-to-onc correspondence can be established between link 
variables and plaquettes. A convenient way of exploiting this fact consists in 
fixing the gauge [94] 

U,,o = 1 (4.1) 

(the lattice version of the temporal gauge Aq = 0). An extremely important 
consequence of the gauge choice (4.1) emerges from considering the gauge-fixed 
form of the single-plaquettc contribution to the lattice action: 

Tr (C/.,oC/.+o,iC/i+i,oC^ii) ^ U,+o,iUli. (4.2) 

This is nothing but the single-link contribution to the one-dimensional lattice 
action of a principal chiral model whose links lie along the direction. When 
considering invariant expectation values (Wilson loops), we then recognize 
that they can be reduced to contracted products of tensor correlations of 
variables defined on decoupled one-dimensional models. As a consequence, 
YM2 factorizes completely into a product of independent chiral models labeled 
by their 1 coordinate. Not only the partition function, but also all invariant 
correlations can be systematically mapped into those of the corresponding 
chiral models. The area law for non self-interacting Wilson loops in YM2 and 
the exponential decay of the two-point correlations in one-dimensional chiral 
models are trivial corollaries of these results [94] . 

The above considerations allow us to focus on the prototype model defined by 
the action 

S = -N^Tr{U,Ul, + UlU,^,), (4.3) 

i 

where i is the site label of the one-dimensional lattice. By straightforward 
manipulations we may show that the most general nontrivial correlation one 
really needs to compute involves product of invariant operators of the form 

Tr(t/ot//)^ (4.4) 

where / plays the role of the space distance, and A; is a sort of "winding 
number" . 
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An almost trivial corollary of the above analysis is the observation that YM2 
and principal chiral models in one dimension enjoy a property of "geometriza- 
tion", i.e., the only variables that can turn out to be relevant for the complete 
determination of expectation values are the single-plaquette (single-link) av- 
erages of products of powers of moments [95] 

ulT^mirp (4.5) 



and the geometrical features of the correlations (in YM2, areas of Wilson loops 
and subloops; in chiral models, distances of correlated points), such that all 
coupling dependence is incorporated in the expectation values of the quantities 
(4.5). This result is sufficiently general to apply not only to the Wilson action 
formulation, but also to all "local" actions such that the interaction depends 
only on invariant functions of the single-plaquette (single-link) variable, i.e., 
any linear combination of the expressions appearing in Eqs. (2.8), (2.9) [95- 
97]. 

In order to proceed to the actual computation, it is convenient to perform a 
change of variables, allowed by the invariance of the Haar measure, parame- 
terizing the fields by 

Vi = Ui^.Ul; (4.6) 

the action (4.3) explicitly factorizes into 

S^-Nj2T^r{Vi + Vi^). (4.7) 
I 



It is now easy to get convinced that in the most general case a Wilson loop 
expectation value (correlation function) can be represented as a finite product 
of invariant tensors, each of which is originated by a single-link integration of 
the form 



JdVif{Vi)exp 


'N(3Tr{Vi + Vi^)] 


J dVi exp 


'NpTr{Vi + Vi^)] 



(4.8) 



where /(VJ) is any (tensor) product of VJ's and s, and the only nontrivial 
contributions to the full expectation value come from integrations extended to 
plaquettes belonging to the area enclosed by the loop itself (in chiral models, 
links comprised between the extremal points of the space correlation) . 
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For sake of definiteness, we may focus on the correlators [98] 



(4.9) 



and find that 



fdVi...dVi{l/N) Tr(\/i...V^)^exp 




Tr{V, + V^) 




'NpTr{Vi + V^) 





(4.10) 



This problem can be formally solved for arbitrary by a character expansion, 
which we shall discuss in Subs. 4.8. It is however immediate to recognize that 
we are ultimately led to computing the general class of group integrals whose 
form is 



y dy n (Ti- V'^y exp [N(3 Tr{V + V^) 



(4.11) 



(where the product runs over positive and negative values of k), and in turn 
it is in principle an exercise based on the exploitation of the result for the 
external fiend single-link integral introduced in Eq. (3.1). 

By the way, integrals of the form (4.11) can easily be expressed as linear 
combinations of integrals belonging to the class 



/ dVx(x){V) exp [NpTr{V + 0)] , (4.12) 



where A labels properly chosen representations of U(iV). Eq. (4.12) is in turn 
related to the definition of the character coefficients in the character expansion 
of exp[NpTr{V + V^)]. For arbitrary N, as a, matter of principle, X(a)(T^) 
has a representation in terms of the eigenvalues <i).j of the matrix V, while 
Tr{V + l/t) = 2 cos 0, and the measure itself can in this case be expressed 
in terms of the eigenvalues as 

dfx{V)r.l[dcf>iA'{4>,,...,cf>N), (4.13) 



where 



A(0i, ...,0Ar) = det exp ||i(i0j)||, 

A^(0i,...,0^) = n4sin^^^. (4.14) 
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As a consequence, it is always possible to express all U(iV) integrals in the 
class (4.12) in terms of linear combinations of products of modified Bessel 
functions 4(2iV/5), with k < N. 

Let us now come to the specific issue of evaluating the relevant physical quan- 
tities in the large-A^ limit of U(A^) models, and comparing the procedures 
corresponding to different possible approaches. Basic to most subsequent de- 
velopments is the observation that the large- factorization property allows 
us to focus on a very restricted class of interesting correlations, which we label 

by 

Wk=(^^TrV''^ = Wi,k. (4.15) 

The first explicit solution to the problem of evaluating in the large- A^ limit 
was offered by Gross and Witten [94]. To this purpose, they introduced the 
eigenvalue density 

Pi<P) = ^T.^i<f^-<f>^)^ (4-16) 
i 

and considered the group integral defining the partition function of the single- 
link model 

~ / nd0iA'(0i,...,0jv)exp(|2A^/3^cos0,j . (4.17) 



The integral (4.17) can be evaluated in the N ^ oo limit by a saddle-point 
technique [99] applied to the effective action 

2/3 y" p(0) cos(0) d(j)+ J p{<p) p{<i)') log sin d0 d0', (4. 18) 



with the constraint / p(0) d0 = 1. The support of the function p(0) is dynam- 
ically determined. The saddle-point integral equation is 

2/3 sin 0= j d0'p(0')cot^^, (4.19) 



and it is possible to identify two distinct solutions, corresponding to weak and 
strong coupling. When (5 is small, it is easy to find out that 

p(0) = J-(l + 2/3cos0), -7r<0<7r; (4.20) 
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p{(f)) is positive definite whenever /9 < |. When /3 is large, (t)c<T^ and 



p(0) = -cos-^--sm - 



sm 



2 V^c 



2 2/3 



(4.21) 



submitted to the condition (3 > \. Therefore it is possible to identify the 
location of the third-order phase transition [94] : 



(4.22) 



By direct substitution, one finds the values of the free and internal energy (per 
unit link or unit plaquette): 



\2/3-|log2/3-|, /3>|, 



Wi 



1 d F 



2 ' 



(4.23) 



(4.24) 



More generally, one may evaluate Wh from p(0), thanks to the relationship 



Vc 

Wk— J d(j) COS k(f) p{(f)) 



0, 



/3<i, k>2, 



1 



1 



2(5) k-i^'^'-^y /?)' 



(4.25) 



where P^'^^ are the Jacobi polynomials. All Wk are differentiable once in 
/3 — /3c, but their second derivatives are discontinuous. Let us notice that Eqs. 
(4.22), (4.23), and (4.24) are an immediate consequence of Eqs. (3.12) and 
(3.15) for the special choice 



(4.26) 



4-2 The Schwinger-Dyson equations of the two-dimensional Yang-Mills the- 
ory 



It is interesting to obtain the above results from the algebraic approach to 
the Schwinger-Dyson equations of the model. We can restrict Eqs. (2.15) to 
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the set of Wilson loops Ck consisting of k turns around a single plaquette, in 
which case by definition W{Ck) = w^- Formally, the Schwinger- Dyson equa- 
tions do not close on this set of expectation values; however, one may check 
by inspection, using the factorization property of two-dimensional functional 
integral for the Yang-Mills theory, that contributions from other Wilson loops 
cancel in the equations for Wk (this is strictly a two-dimensional property) . As 
a consequence, we obtain the large- A?" relationships [100] 

n 

l3{Wn-l - Wn+l) = WkWn-k , (4.27) 
k=l 

with a boundary condition wq — 1- The solution is found by defining a gener- 
ating function 

oo 

$(t) = ^Wfei*^ (4.28) 

k=0 

and noticing that Eq. (4.27) corresponds to 

m'^-(^-l-wit)^l(^^-^), (4.29) 



which is solved by 



(4.30) 

The condition \wk\ < 1 implies that $(t) is holomorphic within the unitary 
circle. On the boundary of the analyticity domain, t — e^'^ and 

1 

- f [Re$(0) - i] cosk(f)d(f), (4.31) 



and as a consequence we may identify 

Re$(0)-| = p(0). (4.32) 

The positivity condition on p(0) leads to a complete determination of the 
solution, implying either 

<0i=/3, wt = (k>2), (4.33) 
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or p(7r) = 0, which in turn leads to 




-0c < > 0, 




2 ' 



(4.34) 



and (pc is given by Eq. (4.21). It is immediate to check that the resulting 
eigenvalue densities are the same as Eqs. (4.20) and (4.21). 

Let us mention that these methods may in principle be applied to more general 
formulation of the theory based on "local" actions, and in particular Wilson 
loop expectation values can be computed for the fixed-point version of the 
model, corresponding to the continuum action [95]. The fixed-point action in 
YM2 in turn is nothing but the "heat kernel" action [101], discussed in the 
large- A/^ context in Ref. [102]. Large- continuum YM2 is slightly beyond the 
purpose of the present review. We must however mention that in recent years 
a number of interesting results have appeared in a string theory context. It is 
worth quoting Refs. [103-105] and references therein. 

While the problem of evaluating the more general expectation values Wi_k 
is solved in principle, in practice it is not always simple to obtain compact 
closed-form expressions whose general features can be easily understood. In 
the strong-coupling regime /5 < |, it is not too difficult to determine from 
finite- A?" results the large- limit in the form [98] 



and one may show that the corresponding Schwinger-Dyson equations close 
on the set Wi^k for any fixed I and are solved by Eq. (4.35). As a matter of 
fact, by defining 



one may show that the strong-coupling Schwinger-Dyson equations reduce to 




(4.35) 



00 



(4.36) 



fc=0 



i-i 



(4.37) 



For the interesting values / = 1 and / = 2, Eq. (4.35) reduces to 



$i(t) = l+/3i. 



(4.38) 
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consistent with the strong-couphng solution (4.30), and 

$2(^) = i(0H^4^+l), 



(4.39) 



related to the generating function for the moments of the energy density 



— (Tr- 
1 



/3t(K + v; 



k=l 



(4.40) 



Eq. (4.40) is related to a different approach for solving large- unitary-matrix 
models, based on an integration of the matrix angular degrees of freedom to 
be performed in strong coupling [61,106]. 

The corresponding weak-coupling problem is definitely more difficult. As far 
as we can see, the Schwinger-Dyson equations close only on a larger set of 
correlation functions, defined by the generating function [107] 



<W = ^Tr 
0<k<l, n > 0, 



{Vkr+'Vk+i...vi 



1 



l-tV,...Vi. 



(4.41) 



such that 



^i{t) = l + tDfUt). 



(4.42) 



The explicit form of the equations is 



n-l 



j=0 



+p =0, i<k<i. 



)(0 



(4.43) 



When / = 1, 2 it is possible to find explicit weak-coupling solutions, but the 
general case / > 2 has not been solved so far. 

More about the calculability of Wilson loops with arbitrary contour in two- 
dimensional U(oo) lattice gauge theory can be found in Ref. [108]. The corre- 
sponding continuum calculations are presented for arbitrary U(A^) groups in 
Ref. [109]. 



39 



4-3 Large-N properties of the determinant 



It is quite interesting to apply the results of Subs. 3.3, concerning the proper- 
ties of the determinant, to YM2 and principal chiral models in one dimension. 
Exploiting the factorization of the functional integration and the possibility of 
performing the variable change (4.6) in the operators as well as in the action, 
we can easily obtain the relationship 



A, = det 



UoUl] = det [Vi...Vi\ = det Vi... det VJ, 



(4.44) 



and, as a consequence, 

(Ai) ^ {detV^y . (4.45) 



The problem is therefore reduced to that of evaluating (det V) in the single- 
plaquette model. It is immediate to recognize from Eqs. (3.31) and (3.32) 
that 




(4.46) 
(4.47) 

Apparently, this expectation value acts as an order parameter for the phase 
transition between the weak- and strong-coupling phases. More precisely, ac- 
cording to Green and Samuel [110,111], one must identify the order parameter 
with the quantity 

{Aif" (4.48) 



and notice that 



(A;) — s> 1 in weak coupling, (4.49) 

(A;)^^'^ — >• exp(— cr/) in strong coupling, (4.50) 

where a acts as a U(l) "string tension". Eqs. (4.49) and (4.50) generahze to 
higher dimensions, when replacing / with the (large) area of the corresponding 
Wilson loop. Notice that the weak-coupling result is consistent with the de- 
coupling of the U(l) degrees of freedom from the SU(iV) degrees of freedom, 
and with the interpretation of U(l) as a free massless field. 



40 



It is therefore interesting to compute 



<J^-j^log{detV) 



(4.51) 



in the case of the single-matrix model; this requires taking the laxge-N limit 
only after the strong-coupling calculation of (det V) has been performed. Since 
the technique of evaluation of a has some relevance for subsequent develop- 
ments, we shall briefly sketch its essential steps. Standard manipulations of 
the single-hnk integrals for finite N allow to evaluate 



AmAP)- J dVew[N(3Tr{V + V^)] {det ^ det \\h-i-m{'^N(3)\\. 

(4.52) 

These quantities can be shown to satisfy the recurrence relations [112] 

(4.53) 



Willing to compute expectation values, we define 

Am.,N 



A^,iv(/3) = ((detn 



A, 



0,N 



(4.54) 



Eq. (4.53) implies that 



^m,N ~ ^m+l,N^m-l,N — ■^m.W-l Am.AT+l (1 — ■^i^at)- 



(4.55) 



Since all Am,i are known, it is possible to reconstruct all A^^^v from Eq. (4.55) 
once Ai_iv is determined. Now Ai^at is exactly (det V) , and it is possible to 
show that it obeys the following second-order differential equation [113] 



1 d d ^ 1 

■ S — Ai AT + 



s ds " ds ' 1 - 

+ (1 - A^ )A^_^ = 0, 



d \' 



ds 



Ai. 



N 



(4.56) 



where s = 2Nf3. Eq. (4.56) can be analyzed in weak and strong coupling and 
in the large- A?" limit. In particular the weak-coupling 1/N expansion leads to 



Ai. 



N 



1 V 



5/2 



2P N'12SP^\'-W) """l^J- 



1 
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thus confirming Eq. (4.46), while in strong couphng one may show that 

Ai,^ = Jn{2N(5) + 0(/53^+2) J^(2iV/?), (4.58) 



AT— >oo 



where is the standard Bessel function, whose asymptotic behavior is well 
known. As an immediate consequence, we find 

-<7 = yr^-logi±^^i^, (3<\. (4.59) 



This result was first guessed by Green and Samuel [111], and then explicitly 
demonstrated in Ref. [113]. 



4-4 Local symmetry breaking in the large-N limit 



Another interesting application of the external-field single-link integral to the 
large- limit of two-dimensional Yang- Mills theories is the study of the possi- 
bility of breaking a local symmetry, as a consequence of the thermodynamical 
nature of the limit. If we introduce an infinitesimal explicit U(A') symmetry 
breaking term in the action [114] 

S = -(3N \TrV + JNV'^ + h.c.l , (4.60) 



corresponding to replacing 

Air^^P[6im + NJ6ijSmi] (4.61) 

in Eq. (3.1), we find that the eigenvalues of AA^ are 

2^1,2 = [l + ^N'^J^ ± I V J^Ar4 + Ajm^] , 

l>2. (4.62) 



When taking the large- A^ limit of the free energy, we find 

lim — — - 

Ar-»oo 



hm l^ = Fo(/?) + 2/3|J|, (4.63) 



and in the limit J ^ 0^ we then find 

'Rc\/*A = ±1. (4.64) 
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We therefore expect that, for finite k, the U{k) global symmetries of large- 
N chiral models and U(A;) gauge symmetries are broken in any number of 
dimensions [114]. This phenomenon cannot occur for any finite value of N in 
two dimensions. 



4-5 Evaluation of higher- order corrections 



In the context of large- two-dimensional Yang-Mills theory, it is worth men- 
tioning that it is possible to compute systematically higher-order corrections 
to physical quantities in the powers of 1/N'^. It is interesting to notice that 
the weak-coupling corrections to the free energy [115] (see also [82]) 



1 1 1 / 1 ' 

A log log 1 

12 12 ^ 8^1 2/3 , 



+ 




(4.65) 



1024/34 



(4.66) 



where A = 0.24875..., arc well defined, but become singular when /3 ^ |. In 
turn, when evaluating higher-order corrections in the strong-coupling phase, 
one finds out that there are no corrections proportional to powers of 1/iV, while 
there are contributions that fall off exponentially with large A^, as expected 
from the general arguments discussed in Subs. 3.4 in connection with the 
appearance of the DeWit-'t Hooft poles. 

Let us however mention that Eqs. (4.55) and (4.56) are also the starting point 
for a systematic 1/N expansion of the free energy in the weak-coupling regime, 
alternative to Goldschmidt's procedure. The basic ingredient is the observation 
that, defining the free energy at finite N by 



F^(/3)=logAo,^(/3), 



(4.67) 



one may show that 



— (logFjv - logFjv-i) 
ds 




(4.68) 
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and this allows for a systematic reconstruction of Fn, whose strong-coupling 
form is [112] 

oo 

Fr,{(3) = N'/3' - Y: kJl+k^P) + 0(/?^^+^). (4.69) 
jfc=i 

4.6 Mixed- action models for lattice YM2 

Another instance of the problem of the single-link integration for matrix fields 
in the adjoint representation of the full symmetry group occurs in the discus- 
sion of the so-called "mixed action" models. Consider the following single-link 
integral [116], resulting from a different formulation of lattice YM2, 

Z{l3u /9a) = / dt/exp [Npi Tr(C/ + U^) + I3^\ Tr U\^] . (4.70) 

It is possible to show that, in the large- iV limit, the corresponding free energy 
can be obtained by the same saddle-point technique presented in Subs. 4.1, i.e., 
by introducing a spectral density p{9) for the eigenvalues of U . This spectral 
density turns out to be precisely the same as the one obtained when /3a = 0, 
if one simply replaces /3f by an effective coupling 

/3eff = A + /9aWl(/3eff), (4.71) 

where w\ can be evaluated in terms of p{9) as 

wi(/3eff) = / d^cos^p(^). (4.72) 

Eq. (4.72) is a self-consistency condition for wi, which allows a determination 
of /3cfr(/3f, /3a). Finally, by substitution into the effective action, one finds the 
relationship 

F(/3f, /3a) = F(/3efr(/3f, /3a), 0) - P^wUPM, /3a)), (4.73) 

where F(P,0) is nothing but the free energy obtained in Subs. 4.1. 

The strong- and weak-coupling solutions are separated by the line 2/3f-|-/3a — 1. 
In strong coupling one obtains 



/3eff = = 
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while in weak coupling 



A + /?a+V(A + /?a)'-/?a 

/9a - A + V(A + /9a)'-/3a 



I' — Pf -\ lo 



A + /3a+ V(A + /5a)2-/3a 



(4.75) 



It may be interesting to quote explicitly the limiting case — 0, where [79] 



Z(0,/?a)= J dt/exp/3a|TrC/| 



0, 



/5a <1, 



(4.76) 



One may actually show that, in any number of dimensions, a lattice gauge 
theory with mixed action [117-119] (a trivial generalization of Eq. (4.70)) is 
solved in the large- limit in terms of the solution of the corresponding theory 
with pure Wilson action; Eqs. (4.71) and (4.73) hold as they stand, and 

(4.77) 

/3f=/3eff, /9a=0 



More about the large-A^ behavior of variant actions can be found in Refs. 
[120-122]. Different kinds of variant actions have been studied in the large- 
limit in Refs. [123-125]. 



4-7 Double-scaling limit of the single-link integral 



In the Introduction, we mentioned that one of the most interesting phenomena 
related to the large- limit of matrix models is the appearance of the so-called 
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"double-scaling limit" 

^ oo, 



N^^^'{gc-g)^ const, (4.78) 
9 9c, 



where g is a (weak) coupling related to the inverse of j3. We already discussed 
the general physical interpretation of this limit as an alternative description of 
two-dimensional quantum gravity and its relationship to the theory of random 
surfaces. Here we only want to consider the double-scaling limit properties for 
those simple models of unitary matrices that can be reformulated as a single- 
link model (cfr. Ref. [126]). 

This specific subject was pioneered by Periwal and Shevitz [127], who discussed 
the double-scahng limit in models belonging to the class 



J dUexp [NpTrV{U + U^)] , (4.79) 



where V{U) is a polynomial in U. Because of the invariance of the measure, 
Eq. (4.79) can be reduced to 

Zj,r^ J d(/)i|A(e'<^\...,e'^^)|2exp[7V/3EiV(2cos(/),)], (4.80) 



and solved by the method of orthogonal polynomials. One starts by defining 
polynomials 

n-l 

Pn{z) = + E (^k,nz\ (4.81) 
k=0 



that satisfy 

dz 



hn 5mn , (4.82) 



where the integration runs over the unit circle, and moreover obey the recur- 
sion relation 

P^+^{z) = zP^iz) + Kz^Pn (^) , ^ = 1 - (4.83) 

where Rn = ao,n+i- As a corollary, 

Z^^N\^(1-RUY~\ (4.84) 
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and one may show that 



(n + 1) (/in+i -hn)= exp 



X (l-l)p„+i(^)P„Q, (4.85) 



which in turn leads to a nonhnear functional equation for 

The simplest example, corresponding to YM2, amounts to choosing V' — 1, 
obtaining 

(n + l)Rl = N(5Rr,{Rn+i + Rn-i){l - Rl), (4.86) 



and in the large- limit, setting n — N and Rn — R, we obtain the limiting 
form 

R^ ^2(3R^{1- R^), (4.87) 



showing that 13c — \ (degeneracy of solution Rc — ^). One may now look for 
the scahng solution to Eq. (4.86) in the form 

R^-Rc^Rn^ N-'^f [NP{gc -g)], g ^ ^, (4.88) 



where is related to the second derivative of the free energy. This is a con- 
sistent Ansatz when 

11=1, p=l (4.89) 

leading to the equation 

-2xf + 2f = f", x = NP{gc-g). (4.90) 

In the case V' = 1 -\- \u, one finds the equation 

i = -2(l-i?2)(-l- A + 3Ai?2), (4.91) 

which reduces to 1/13 = |(1 — i?^) when A = |. A scaling solution to the 
corresponding difference equation requires /x = | and p = |. When V' = 
1 -I- Xiu + A2W^, multicriticality sets at Ai = — | and ^2 — j^, and = 
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y(1 — R^), leading to the exponents = y and p = |. Rather general results 
can be obtained for an arbitrary order k of the polynomial V: /J, — l/{2k+ 1), 
p = 2k/ {2k + 1), and c = 1 - Q/{k{k + 1)). 

The double-scahng limit can also be studied in the case of the external-field 
single-link integral [128], and it was found that its critical behavior is simple 
enough to be identified with that of the k = 1 unitary-matrix model. In the 
language of quantum gravity, the only effect of introducing A^^ real parameters 
Aij is that of renormalizing the cosmological constant, without changing the 
universality class of the critical point. 

A few interesting features of the double-scaling limit for the k = 1 model are 
worth a more detailed discussion [21]. In particular let us recall that, according 
to Eq. (4.89), 




and therefore 71 = 3, implying c = —2. We may now reinterpret the double- 
scaling limit of matrix models as a finite-size scaling with respect to the "vol- 
ume" parameter in a two-dimensional NxN space. As a consequence, we 
obtain relationships with more conventional critical exponents through the 
identification 71 = 2i/, which in turn by hyperscaling leads to a determination 



3 
2 

a = —1. The result a = —1 can be easily tested on the solution of the model 



of the specific heat exponent a = 2(1 — u). Numerically we obtain = | and 



with Pc — 21 consistent with a negative critical exponent a — —1. 

It is also interesting to find tests for the exponent z/, especially in view of 
the fact that the most direct checks are not possible in absence of a proper 
definition for the relevant correlation length. Numerical studies have been per- 
formed by considering the partition function zero /?o closest to the transition 
point /9c = |, finding that the relationship 

Im/3o oc AT-V'^ (4.94) 



is rather well satisfied even for very low values of A^; at A^ > 5, it is valid 
within one per mille. Another test concerns the location of the peak in the 
specific heat in U(A^) models, whose position Ppea.k{N) should approach /3c 
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with increasing N. Finite-size scaling arguments predict 

/3peak(iV)^/5e + a7V-V-, (4.95) 

and large- A?" results are very well fitted by the choice = |, a = 0.60 [129]. 

4-8 The character expansion and its large-N limit: SU(A^) vs. U(A/") 

The general features of the character expansion for lattice spin and gauge 
models have been extensively discussed by different authors. In particular, 
Ref . [77] , besides offering a general presentation of the issues, presents tables of 
character coefficients for many interesting groups, including U(oo) = SU(oo), 
for the Wilson action. Let us therefore only briefiy recall the fundamental 
points of this approach, which is relevant especially in the analysis of the 
strong-coupling phase and of the phase transition. 

In Sect. 2 we classified the representations and characters of U(A'") groups. Be- 
cause of the orthogonality and completeness relations, every invariant function 
of V can be decomposed in a generalized Fourier series in the characters of V. 
Let us now consider for sake of definiteness chiral models with action given by 
Eq. (2.10); extension to lattice gauge theories is essentially straightforward, 
at least on a formal level. We can replace the Boltzmann factor corresponding 
to each lattice link by its character expansion: 



exp {pNTr [uM+t^ + U.+^.Ul\ } 



(4.96) 



where the sum runs over all the irreducible representations of \](N), F(f3) is 
the free energy of the single-link model 



F{(3) = ^log / dVe^v[N(3TT{V + V^)] = ^ logdet ||/,_,(27V/?)||, 

(4.97) 

and 5(r)(/3) are the character coefficients, defined by orthogonality and repre- 
sentable in terms of single-link integrals as 



d{r)Z[r){P) = {X{r){V) 



det 11/ 



det||J,_,(2iV/5)|| ' 



(4.98) 
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with A defined by Eq. (2.2). We may notice that, for any finite N, Z(^r){P) are 
meromorphic functions of f3, with no poles on the real axis, which is relevant 
to the series analysis. However, singularities may develop, as usual, in the 
large-A^ limit. Eqs. (4.96) and (4.97) become rapidly useless with growing A^. 
However, an extreme simplification occurs in the large- A?" limit, owing to the 
property 



Ti-\-\ Ti—\ '- 



(4.99) 



where n+ — J2ih, ''T'- — and (7(;) is the dimension of the representation 

(l) of the permutation group, which in turn can be computed explicitly as 

d(i,m) can be parameterized by 

d{i,m) = — r — :<7(i)<^{m) Ca m), (4.101) 
n+! n_! 

where C(/^^) can be expressed as a finite product: 



_^ (N-t-i + h)\ ^ {N-s-j+m,)\ 

ti/ii {N + l-i-j)\ ^ ^ 

allowing for a conceptually simple 1/N expansion. These results are comple- 
mented with the result 

F{(3) = /?2 + 0(/?2^+2) (4.103) 



and with the unavoidable large- A?" constraint /3 < |. 
The character expansion now proceeds as follows. 

We notice that, thanks to Eq. (4.99), only a finite number of nontrivial rep- 
resentations contributes to any definite order in the strong-coupling series ex- 
pansion in powers of (3, and each lattice integration variable can appear only 
once for each link where a nontrivial representation in chosen. A systematic 
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treatment leads to a classification of contribTitions in terms of paths (surfaces 
in a gauge theory) along whose non self-interacting sections a particular rep- 
resentation is assigned. Self-intersection points are submitted to constraints 
deriving from the orthogonality of representations and their composition rules. 

In the case of chiral models, all relevant assignments can be generated by 
considering the class of the lattice random paths satisfying a non-backtracking 
condition [130]. 

Once all nontrivial configurations are classified and counted, one is left with 
the task of computing the corresponding group integrals. Only integrations at 
intersection points are nontrivial, since other integrations follow immediately 
from the orthogonality relationships. Unfortunately, no special computational 
simplifications occur in the large- A?" limit of group integrals. 

Apparently, the character expansion is the most efficient way of computing the 
strong-coupling expansion of lattice models. In particular, very long strong- 
coupling series have been obtained in the large-A^ limit for the free energy, 
the mass gap, and the two-point Green's functions of chiral models in two 
and three dimensions (for the free energy, 18 orders on the square lattice, 26 
orders on the honeycomb lattice, and 16 orders on the cubic lattice; for the 
Green's functions, 15 orders on the square lattice, 20 orders on the honeycomb 
lattice, and 14 orders on the cubic lattice). The analysis of these series will be 
discussed in Sect. 7. 

Before leaving the present subsection, we must make a few comments con- 
cerning the relationship between SU(A^) and U(A^) groups. We already made 
the observation that when iV — >■ oo there is essentially no difference between 
SU(A^) and U(A'') models, at least when considering operators not involving 
the determinant. In order to explore this relationship more carefully, we may 
start as usual from the expression of the single-link integral (3.1). 

Representations of Z{A^A) in the SU(A^) case can be obtained [131] in terms 
of the eigenvalues Xi of A^A and of 9, defined in Eq. (3.23). Introducing the 
Vandemonde determinant 

A(Ai,...,A^) = n(Ai - AO =det||A}-i, (4.104) 

j>i 

one obtains 

z{A^A) - ^ ( n' I) / n ^ (e + Ne) 
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X 



A{2y/x^, ...,2y/x]^) A(cos 01, ...,cos 0jv) 



exp 



k 

(4.105) 



or alternatively 



AT-l 



X ^- exp 

A(a;i, ...,xn) 



k 



(4.106) 



The only difference between SU(A^) and U(A^) is due to the presence of the 
(periodic) delta function S {J2i 4>i + -^^)) introducing the dependence on 9 cor- 
responding to the constraint detU = 1. A formal solution is obtained by 
expanding in powers of e^^^: 



Z{A^A)= J2 e'^™^det||;.rVi-H(2^0ll 



.k=l 



A(^?,...,4)' 

(4.107) 



where Zi = ^/Xi. Eq. (4.107) in turn leads to the following representation of 
the free energy for the SU(A^) single-link model: 



Fjv(/3,e) = log £ Am,N(P)e 



iNme 



(4.108) 



where for convenience we have redefined the coupling: [3 — > (3e^^. Eq. (4.108) 
is useful for a large- mean-field study [112], but it is certainly inconvenient 
at small N, where more specific integration techniques may be applied. 

We mention that a large- A^ analysis of Eq. (4.108) for ^ = leads to 



Fjv(/3, 0) = N^P"" + 2M2N(3) - 2 Jjv-i(2Ar/3) Jjv+i(2Ar/3) 

oo 

-E^4+.(2Ar/5)-hO(/?^^). 



(4.109) 



k=l 



It is also possible to establish a relationship between SU(A^) and U(A") groups 
at the level of character coefficients. Thanks to the basic relationships 



(4.110) 
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holding in U(A'^), one may impose the condition detU = 1 in the integral 
representation of the character coefficients and obtain 




oo 



oo 



z{r, s) 
z{0,sy 



(4.111) 



where, by definition, for U(A'") groups 



z{0, s) = (det , z{r, s) = (det U'X{r){U)) . (4.112) 



These relationships arc the starting point for a systematic implementation of 
the corrections due to the SU(A^) condition in the 1/A^ expansion of U(A^) 
models [52,132]. A peculiarity of the SU(A^) condition can be observed in the 
finite- behavior of the eigenvalue density function p{(f>,N), which shows a 
non-monotonic dependence on (p, characterized by the presence of N peaks. 
This is already apparent in the /3 ^ limit of the single-link integral, where 



5 Chiral chain models and gauge theories on polyhedra 

5. 1 Introduction 

The use of the steepest-descent techniques allows to extend the number of the 
unitary-matrix models solved in the large- limit to some few unitary-matrix 
systems. The interest for few-matrix models may arise for various reasons. 
Their large- A^ solutions may represent non-trivial benchmarks for new meth- 
ods meant to investigate the large- A^ limit of more complex matrix models, 
such as QCD. Every matrix system may have a role in the context of two- 
dimensional quantum gravity; indeed, via the double scaling limit, its critical 
behavior is connected to two-dimensional models of matter coupled to gravity. 
Furthermore, every unitary-matrix model can be reinterpreted as the gener- 
ating functional of a class of integrals over unitary groups, whose knowledge 
would be very useful for the strong-coupling expansion of many interesting 
models. 

This section is dedicated to a class of finite-lattice chiral models termed chain 



[129] 



(4.113) 
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models and defined by the partition function 

. L r ^ 

Zl^JU dUi exp 5^ IV {UiUl, + U}u,+^ 



i=l 



i=l 



(5.1) 



where periodic boundary conditions are imposed: = Ui. 

Chiral chain models have interesting connections with gauge models. Fixing 
the gauge Aq = 0, YM2 on a x L lattice (with free boundary conditions in 
the direction of size K) becomes equivalent to K decoupled chiral chains of 
length L. 

Chiral chains with periodic boundary conditions enjoy another interesting 
equivalence with lattice gauge theories defined on the surface of polyhedra, 
where a link variable is assigned to each edge and a plaquette to each face. 
By choosing an appropriate gauge, lattice gauge theories on regular polyhedra 
like tetrahedron, cube, octahedron, etc., are equivalent respectively to periodic 
chiral chains with L = 4,6, 8, etc. [70]. 

The thermodynamic properties of chiral chains can be derived by evaluating 
their partition functions. Free-energy density, internal energy, and specific heat 
are given respectively by 

FL^j^logZ,, (6.2) 



When L ^ 00, Zl can be reduced to the partition function of the Gross- 
Witten single- link model, and therefore shares the same thermodynamic prop- 
erties. In particular, the free energy density at = 00 is piecewise analytic 
with a third-order transition at (3^ = \ between the strong-coupling and weak- 
coupling domains. Furthermore, the behavior of Coo around {3^ can be charac- 
terized by a specific heat critical exponent a = —1. It is easy to see that the 
L — 2 chiral chain is also equivalent to the Gross- Witten model, but with P 
replaced by 2/3; therefore Pc — ^ and the critical properties are the same, e.g., 
a — —1. 
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5.2 Saddle-point equation for chiral L-chains 



The strategy used in Refs. [70,133] to compute the N — oo solutions for chiral 
chains with L < 4 begins with group integrations in the partition function 
(5.1), with the help of the single-link integral, for all Ui except two. This leads 
to a representation for Zl in the form 



(5.5) 



suitable for a large-A^ steep est- descent analysis. Since the integral depends 
only on the combination UV\ changing variable to 6j, e^^J being the eigenval- 
ues of UV'^ , leads to 

Zl^ J nd^i|A(^i,...,Mrexp [N'sii\9k)] (5.6) 



where -tt < 9j < tt, A{9i, ...,e^:) = det||Ajfe||, A^-^ = e'^'^^ In the large- iV 
limit, Zl is determined by its stationary configuration, and the distribution of 
9j is specified by a density function Pl{9), which is the solution of the equation 



/ d<p pl{<P) cot ^ + j^Sii^ {9, pl) = 0, 



(5.7) 



with the normahzation condition 



TT 

/ Pl{0) 



d9^1. 



(5.8) 



For L = 2, ^2 is already in the desired form with 

Si^ = 2/3^ Tr (U,UI + UlU^) , 



(5.9) 



and the large- A?" eigenvalue density P2{9) of the matrix UiUl satisfies the Gross- 
Witten equation 



j d(/) p2{(t>) cot - A(3 sin ^ = 0, 



(5.10) 



which differs from that of the infinite-chain model only in replacing /5 by 2/3. 
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5.3 The large-N limit of the three-link chiral chain 



In the L — 3 chain model, setting U — Ui and V — U2, S^g is given by 



exp 



= exp [2N(3ReTrUV^] J dC/gexp [2N(3ReTr AUl] , (5.11) 



where A — U -\-V. Recognizing in the r.h.s. of (5.11) a single-hnk integral, one 
can deduce that the large- limit of the spectral density P3{9) of the matrix 
UV^ satisfies the equation 



2/5 (sin ^ + sin i^) - j ^(p pz{<i)) 



d-(t) 1 
cot — z h 



sin \9 



2 2 cos \e + cos \(t) 



0, 

(5.12) 



with the normalization condition / p^{9) d6' = 1. In order to find a solution 
for the above equation, one must distinguish between strong-coupling and 
weak-coupling regions. 



In the weak-coupling region the solution of Eq. (5.12) is 



P^{0) = -cos- 

TT 4 



2cos- Wl 

2 V 3/3 



2cos- - 24/1 - ^ 
2 V 3/3 



1/2 



(5.13) 



for 



1^1 < 9r — 2arccos wl 



3(3 



(5.14) 



and ^3(6*) = for < 16*1 < n. This solution is valid for /3 > /3c = |, indicating 
that a critical point exists at /3c = |. Similarly one can calculate P3{9) in the 
strong-coupling domain /3 < /3c [70,133,134] finding: 



Ps{0)-^{y{e) + i 2 , 



X 



(y{0) + x/S) (y{e) + 74+^) 



1/2 



(5.15) 



where 



y{e) = W4cos2- + c, 



(5.16) 
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and the parameter c is related to (3 by the equation 

l + VS+|c+(l-iVc)\/4T^=^. (5.17) 

At P — Pc, c — and therefore 

P3(^)crit = ^ ( 2C0S- j COS-, (5.18) 

in agreement with the critical limit of the weak-coupling solution (5.13). 

Since P3(7r) > for /3 < /3c and psi^n) = for P > [3c, the critical point Pc 
can be also seen as the compactification point for the spectral density P3{0), 
similarly to what is observed in the Gross- Witten model. 

5.4 The large-N limit of the four-link chiral chain 



For -L = 4, setting U — Ui and V — C/3, S^g is given by 



exp (N^SiS) = J d^/aexp (2iV/3 Re Tr 

X J dC/4exp (2NpReTrAUl) , (5.19) 

where again A = U -\-V. The large- limit of the spectral density ^4(6*) of the 
matrix UV^ must be solution of the equation 



4/3 sin jd(t)pi{(t)) 



e-d) sin \e 

cot ' ^ 



2 cos \9 + cos 1^ 



= 0, (5.20) 



satisfying the normalization condition / p4(^)d^ = 1. 



In order to solve Eq. (5.20) one must again separate weak- and strong-coupling 
domains. In the weak-coupling region the solution is 



p^{e) = — Wsin^ ^ - sin^ - for < ^ < < tt. 



TT 



2 2 - ^-c^ (5 21) 



P4(^) =0 for ^c < ^ < TT, 



with implicitly determined by the normalization condition /^^^ p4(^)d^ = 1. 
The solution (5.21) is vahd for P > Pc — l^r, since the normalization condition 
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can be satisfied only in this region. gVr is then a point of non-analyticity 
representing the critical point for the transition from the weak to the strong- 
coupling domain. 



In the strong-coupling domain P < Pc — h'^ one finds 




(5.22) 



where A is determined by the normalization condition J^^p4^{6)d6 = 1. The 
strong- and weak-coupling expressions of p4{9) coincide at (3c' 



Notice that again the critical point f3c = l^r represents the compactification 
point of the spectral density P4{9); indeed P4(7r) > for /9 < P^, and P4(7r) = 



5. 5 Critical properties of chiral chain models with L < 4 

In the following we derive the N = oo critical behavior of the specific heat in 
the models with L = 3,4, using the exact results of Subs. 5.3 and 5.4. 

From the spectral density pz{0), the internal energy can be easily derived by 
U-s = J d9p3{9)cos9. One finds that U3 is continuous at Pc- In the weak- 
coupling region P > Pc — I, 




(5.23) 



for P > Pc. 




(5.24) 



Close to criticality, i.e., for < /3//3c 



1 < 1, 



^3=g-^(/^-/^c//V0(/^ 



(5.25) 



In the strong-coupling region, one finds 




72 2V3 



17 1 



iPc - pT + 0{Pc 



(5.26) 
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for < 1— /9//9c ^ 1- Then the weak- and strong-coupUng expressions of 
C3 show that the critical point P^. = | is of the third order, and the critical 
exponent associated with the specific heat is a = — |. 

In the L — 4: case, recalling that p4:{9) is the spectral distribution of C/1C/3, one 
writes 



r 6 r [ 6 

8/3 1 d6 pi{9) cos - - / d9 dcj) pa{9) pa{4>) log f cos - + cos ■ 

- ^-log2/3+ -^d^d(/)p4Wp4(0)logsin2^ 



(5.27) 



Observing that, since Pi{9) is a solution of the variational equation SF^/Sp^ 
0, the following relation holds 



dF4 _ 9F4 

d^ ~ w 



(5.28) 



one can easily find that 



C/4 = -^ + / d^p4(^)cos 



(5.29) 



In this case, the study of the critical behavior around (3c = |7r is slightly 
subtler, since it requires the expansion of elliptic integrals F{k) and E{k) 
around k = 1. Approaching criticality from the weak-coupling region, i.e., 
when /3 ^ one obtains 



where 5^ ~ /3 — /9c, apart from logarithms. For /3 — > /3j 



where 6^ j3c — apart from logarithms. A comparison of Eqs. (5.30) and 
(5.31) leads to the conclusion that the phase transition is again of the third 
order, with a specific heat critical exponent a — 0~ . 

In conclusion we have seen that chain models with L = 2, 3, 4, 00 have a third- 
order phase transition at increasing values of the critical coupling, /3c = ^, |, 
Itt, I respectively, with specific heat critical exponents a — —1, — |, 0~, — 1 
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respectively. It is worth noticing that a increases when L goes from 2 to 4, 
reaching the hmit of a third order critical behavior, but in the large-L limit it 
returns to a = — 1. 

The critical exponent describing the double-scaling behavior for A?" — > oo 
and j3 ^ j3c, can then be determined by the two-dimensional hyperscaling 
relationship 2z/ = 2 — a. This relation has been proved to hold for the Gross- 
Witten problem, and therefore for the L = 2 and L = oo chain models, where 
it is related to the equivalence of the corresponding double scaling limit with 
the continuum limit of a two-dimensional gravity model with central charge 
c = — 2. It is then expected to hold in general for all values of L. At L = 4, the 
value V = 1 has been numerically verified, within a few per cent of uncertainty, 
by studying the scaling of the specific heat peak position at finite N . Notice 
that the exponents a — Q~ , v — 1 found for L = 4 correspond to a central 
charge c — 1. 

5.6 Strong- coupling expansion of chiral chain models 

Strong-coupling series of the free energy density of chiral chain models can be 
generated by means of the character expansion, which leads to the result 

Fl{(5) = F{(5) + Fl{(5), (5.32) 

where F[j3) is the free energy of the single unitary-matrix model, 

^L = -^logE4)4)' (5-33) 

(r) 

denotes the sum over all irreducible representations of U(A^), and d^r) 
and Z{r){P) are the corresponding dimensions and character coefficients. The 
calculation of the strong-coupling series of FlI/S) is considerably simplified in 
the large- limit, due to the relationships (4.103) and 

^(,)(/5) = ^(,)/3" + 0(/52^), (5.34) 

where Z(^r) is independent of (3 and n is the order of the representation (r). 
Explicit expressions for d(^r) and Z(^r) were reported in Subs. 4.8. The large-A'^ 
strong-coupling expansion of Fl{i3) is actually a series in f3^, i.e., 

Fi = ^c(n,L)/?"^. (5.35) 

n 
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It is important to recall that the large- iV character coefficients have jumps 
and singularities at /? = | [52], and therefore the relevant region for a strong- 
coupling character expansion is /3 < |. 



Another interesting aspect of the large- limit of chain models, studied by 
Green and Samuel using the strong-coupling character expansion [75], concerns 
the determinant channel, which should provide an order parameter for the 
phase transition. The quantity 



is non-zero in the strong-coupling domain and zero in weak coupling at N — 
oo. Pc may then be evaluated by determining where the strong-coupling eval- 
uation of the order parameter a vanishes. Like the free-energy, a is calculable 
via a character expansion. Indeed 



Green and Samuel evaluated a few orders of the above character expansion, ob- 
taining estimates of (3c from the vanishing point of a. Such estimates compare 
well with the exact results for L = 3,4. In the cases where /3c is unknown, they 
found j3c — 0.44 for L = 5, — 0.47 for L = 6, etc., with f3c monotonically 
approaching the value | with increasing L. 

In order to study the critical behavior of chain models for L > 5, one can also 
analyze the corresponding strong-coupling series of the free energy (5.32) [135]. 
An integral approximant analysis of the strong-coupling series of the specific 
heat led to the estimates Pc — 0.438 for L = 5 and Pc — 0.474 for L = 6, with 
small negative a, which could mimic an exponent a — 0~. For L > 7 a such 
strong-coupling analysis would lead to /3c larger than |, that is out of the region 
where a strong-coupling analysis can be predictive. Therefore something else 
must occur earlier, breaking the validity of the strong-coupling expansion. An 
example of this phenomenon is found in the Gross- Witten single-link model 
(recovered when L — > oo), where the strong-coupling expansion of the N — oo 
free energy is just F{(3) = /3^, an analytical function without any singularity; 
therefore, in this model, f3c — ^ cannot be determined from a strong-coupling 
analysis of the free energy. 

Prom such analysis one may hint at the following possible scenario: as for L < 

4, for L = 5,6, that is when the estimate of ,(3,. coming from the above strong- 
coupling analysis is smaller than | and therefore acceptable. The term F{(3) 
in Eq. (5.32) should be the one relevant for the critical properties, determining 
the critical points and giving a —1 (maybe a = 0~ as in the L — A case). 




(5.36) 



{detUiUl,) = 



J2{r) d{r)z\^r) ^C?(r ,-1)^,-1) 
E(r) '^(r)^M 



(5.37) 
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For L >7 the critical point need not be a singular point of the free energy in 
strong or weak coupling, but just the point where weak-coupling and strong- 
coupling curves meet each other. This would cause a softer phase transition 
with q; = —1, as for the Gross- Witten single-link problem. We expect /3c < | 
also for L > 7. This scenario is consistent with the results of the analysis of 
the character expansion of a, defined in Eq. (5.36). 



6 Simplicial chiral models 



6. 1 Definition of the models 



Another interesting class of finite-lattice chiral models is obtained by consider- 
ing the possibility that each of a finite number of unitary matrices may interact 
in a fully symmetric way with all other matrices, while preserving global chiral 
invariance; the resulting systems can be described as chiral models on (rf — 1)- 
dimensional simplexes, and thus termed "simplicial chiral models" [135,136]. 

The partition function for such a system is: 



exp 



1=1 



1=1 j=i+l 



(6.1) 



Eq. (6.1) encompasses as special cases a number of models that we have already 
introduced and solved; in particular, the chiral chains with L < 3 correspond 
to the simplicial chiral models with d <3. 

One of the most attractive features of these models is their relationship with 
higher-dimensional systems, with which they share the possibility of high coor- 
dination numbers. This relationship becomes exact in the large-ii limit, where 
mean-field results are exact. 

In the large- limit and for arbitrary d a saddle-point equation can be derived, 
whose solution allows the evaluation of the large- free energy 

^<i=4log^<i (6.2) 



and of related thermodynamical quantities. 
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6.2 Saddle-point equation for simplicial chiral models 



The strategy for the determination of the large- A?" saddle-point equation is 
based on the introduction of a single auxiliary variable A (a complex matrix), 
allowing for the decoupling of the unitary matrix interaction: 



^0 



(6.3) 



where 



/a 
HdUidA exp 
i=l 



-Np Tr AA^ + N(3TrAj2U^i 

i 

+N(3 TrA+^t/i- Ar2/5c^ . 



(6.4) 



We are now back to the single- link problem and, since we have solved it in Sect. 
3 in terms of the function W, whose large- limit is expressed by Eq. (3.15), 
we obtain 



Zd^ JdA exp [-N(3TrAA^ + NdW{(3^AA^) - N^(3d 



(6.5) 



It is now convenient to express the result in terms of the eigenvalues Xi of the 
Hermitian semipositive-definite matrix 4(3 A A\ obtaining 



= j diJ,{xi) 



exp 



(6.6) 



The angular integration can be performed, leading to 

dn{xi) = Yl dxi Ylixi - Xjf. 

i i>j 



(6.7) 



The saddle-point equation is therefore 



+ Xi 1 ^ (4: — d)^/r+^i + d^yr+Wj 



2/3 



(6.8) 
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subject to the constraint (needed to define r) 

1 



r = 



1 (strong coupling); 
(weak coupling). 



(6.9) 



The energy 



1 dF. 



2 dp 



(6.10) 



is easily expressed in terms of the eigenvalues: 



(6.11) 



In the large- limit, after a change of variables to Zi = + Xj, we intro- 
duce as usual an eigenvalue density function p{z)^ and turn Eq. (6.8) into the 
integral equation 



2/9 



-d 



j-dz' p{z') 



d-2 



z — z' z -\- z' 



(6.12) 



subject to the constraints 



/ P{z') dz' 



(6.13) 



and 



/ 



/ /X dz' 
p{z') — < 1, 



(6.14) 



with equality holding in strong coupling, where a — \fr. The easiest way of 
evaluating the free energy is the integration of the large- version of Eq. 
(6.11) with respect to /3. 

Very simple solutions are obtained for a few special values of d. When d = 0, 
the problem reduces to a Gaussian integration, and one easily finds that Eq. 
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(6.12) is solved by 



z 



— a 



(6.15) 



47r/3 



and = exp(A^' 



^ log /5) , independent of a as expected. 



When d = 2 we obtain 



(6.16) 




r{P) = (1 - 2/3)^ /3 < i, (6.17) 



and these results are consistent with the reinterpretation of the model as a 
Gross- Witten one-plaquette system. Notice however that the matrix whose 
eigenvalue distribution has been evaluated is not the original unitary matrix, 
and corresponds to a different choice of physical degrees of freedom. This is the 
reason why, while knowing the solution for the free energy of the d = 1 system 
(trivial, non-interacting) and of the d = 3 system (three-link chiral chain), 
we cannot find easily explicit analytic forms for the corresponding eigenvalue 
densities. 

The saddle-point equation (6.12) has been the subject of much study in recent 
times, because it is related to many different physical problems in the context 
of double-scaling limit investigations. In particular, in the range of values 
< d < 4, the same equation describes the behavior of 0(n) spin models 
on random surfaces in the range — 2 < n < 2, with the very simple mapping 
n = d — 2 [137]. In this range, the equation has been solved analytically in 
Refs. [138] and especially [139] in terms of ^-functions. 



6.3 The large-N d = A simplicial chiral model 

The chiral model on a tetrahedron is the first example within the family of 
simplicial chiral models which turns out to be really different from all the 
systems discussed in the previous sections. Exphcit solutions were found for 
both the weak and the strong coupling phases, and they are best expressed in 
terms of a rescaled variable 




(6.18) 
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and of a dynamically determined parameter 



k 



'1- 
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(6.19) 



The resulting expressions, after defining /3p(C) d( = p{z) dz, are 

r'2 _ r2 



PwiO = 



E{ky 



K{k) - V^2-c^yw^^(c^ k) 



(6.20) 



and 



p.(C) 



[E{k) - (1 - A;2)ir(A;)]^ 



X 



-Kik)-Jk^-CJi-CU{C,k) 



where E and 11 are the standard elliptic integrals, and < C < ^- 



(6.21) 



The complete solution is obtained by enforcing the normalization condition, 
which leads to a relationship between (3 and k, best expressed by the equation 



1 



(6.22) 



Criticality corresponds to the limit A; — > 1, and it is easy to recognize that 
both weak and strong coupling results lead in this limit to Pc— \ and 



/?Pe(C)=Cl0g 



1 + C 
1-C 



(6.23) 



Many interesting features of this model in the region around criticality can be 
studied analytically, and one may recognize that the critical behavior around 
Pc = \ corresponds to a limiting case of a third-order phase transition with 
critical exponent of the specific heat a = 0~. In the double-scaling limit lan- 
guage this would correspond to a model with central charge c — 1 and loga- 
rithmic deviations from scaling. The critical behavior of the specific heat on 
both sides of criticality is described by 



TT 



dp fe'-o 36 121og(4/A;') 



+ 



Jog k'^ 



(6.24) 



where k' = y/1 — k'^. 
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6-4 The large-d limit 



By introducing a function defined by 

b / f\ 

f{z)=f^dz', f{z)^-^-, (6.25) 

J Z — Z »oo Z 

a 



analytic in the complex z plane with the exception of a cut on the positive 
real axis in the interval [a,b], we can turn the saddle-point equation (6.12) 
into the functional equation 

^-d^2Ref{z) + {d-2)f{-z). (6.26) 



This equation can be the starting point of a systematic 1/d expansion, on 
whose details we shall not belabor, especially because its convergence for small 
values of d is very slow. It is however interesting to solve the large-d hmit of 
Eq. (6.26) by the Ansatz 

p{z) = d{z-z), (6.27) 
whose substitution into Eq. (6.25) leads to the solution 

' = ^i'^f^)' (6.28) 
z^l, (5d<l. 



The large-d limit predicts the location of the critical point j3c = 1/d, and shows 
complete equivalence with the mean-field solution of infinite- volume principal 
chiral models on a (i/2-dimensional hypercubic lattice. The large-d prediction 
for the nature of criticality is that of a first-order phase transition, with 




6.5 The large-N criticality of simplicial models 



The connection with the double-scaling limit problem naturally leads to the 
study of the finite-/? critical behavior. In the regime < d < 4 one is helped 
by the equivalence with the solved problem of 0(n) spin models on a random 
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surface, which allows not only a determination of the critical value (found 
to satisfy the relationship Pcd = 1), but also an evaluation of the eigenvalue 
distribution at criticality [137]: 

2 7:9 sinh 9u 

Pc{z) = — cos — , 6.30 

nU 2 cosh u 



and 



2 71$ 

0, = ^ tan Y , 



(6.31) 



where 9 and u are defined by the parametrizations 



4 cos^ — = d = — , cosh u = — . (6.32) 
z Pq z 



Unfortunately, the technique that was adopted in order to find the above so- 
lution does not apply to the regime c? > 4, in which case one cannot choose 
Oc = 0. The saddle-point equation at criticality can however be solved numer- 
ically with very high accuracy, and one finds that the relationship 

(3cd = 1 (6.33) 



is satisfied for all d, thus also matching the large-d predictions. The combina- 
tions {qc + bc)/2 and Ocfcc admit a 1/d expansion, and the coefficients of the 
expansion are found numerically to be integer numbers up to order d~^. 

An analysis of criticahty for d > 4 shows that its description is fully consistent 
with the existence of a first-order phase transition, with a discontinuity of the 
internal energy measured by dal/{A{d — 1)), again matching with the large-d 
(mean-field) predictions. 



6.6 The strong- coupling expansion of simplicial models 



There is nothing peculiar in performing the strong-coupling expansion of Eq. 
(6.1). There is however a substantial difference with respect to the case of 
chiral chains discussed in the previous section: because of the topology of 

simplexes, the strong-coupling configurations entering the calculation are no 
longer restricted to simple graphs whose vertices are joined by at most one 
link, and the full complexity of group integration on arbitrary graphs is now 
involved [130]. 
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As a consequence, as far as the simplicial models can be solved by different 
techniques, they may also be used as generating functionals for these more 
involved group integrals, that enter in a essential way in all strong-coupling 
calculations in higher-dimensional standard chiral models and lattice gauge 
theories. 



7 Asymptotically free matrix models 

7.1 Two-dimensional principal chiral models 

Two dimensional SU(A'^) xSU(A'^) principal chiral models, defined by the action 

S^^j d^xTrd^U{x)d,U\x), (7.1) 



are the simplest asymptotically free field theories whose large- A?" limit is a sum 
over planar diagrams, like four dimensional SU(A^) gauge theories. 

Using the existence of an infinite number of conservation laws and Bethe- 
Ansatz methods, the on-shell solution of the SU(A^) x SU(A^) chiral models 
has been proposed in terms of a factorized -S'-matrix [140,141]. The analysis 
of the corresponding bound states leads to the mass spectrum 

M, = M!!H(r!M, l<.<iV-l, (7.2) 

sm(7r/A') 



where is the mass of the r-particle bound state transforming as totally 
antisymmetric tensors of rank r. M = Mi is the mass of the fundamental 
state determining the Euclidean long-distance exponential behavior of the 
two-point Green's function 

G{x)^hTrU{0)U{xy). (7.3) 



The mass-spectrum (7.2) has been verified numerically at A^ = 6 by Monte 
Carlo simulations [142,143]: Monte Carlo data of the mass ratios M2/M and 
M'^/M agree with formula (7.2) within statistical errors of about one per cent. 
Concerning the large- A^ limit of these models, it is important to notice that the 
S*- matrix has a convergent expansion in powers of 1/A^, and becomes trivial, 
i.e., the -S-matrix of free particles, in the large- A^ limit. 
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By using Bcthc-Ansatz techniques, the mass/A-parameter ratio has also been 
computed, and the result is [144] 



M 



V e n/N 



/Stt sin(7r/iV) 



(7.4) 



■MS 



which again enjoys a 1/N expansion with a finite radius of convergence. This 
exact but non-rigorous result has been substantially confirmed by Monte Carlo 
simulations at several values of N [132,145], and its large- limit also by 
N — oo strong-coupling calculations [146,147]. 

While the on-shell physics of principal chiral models has been substantially 
solved, exact results of the off-shell physics are still missing, even in the large- 
A^ limit. When A^ — > oo, principal chiral models should just reproduce a free- 
field theory in disguise. In other words, a local nonlinear mapping should exist 
between the Lagrangian fields U and some Gaussian variables [7]. However, 
the behavior of the two-point Green's function G{x) of the Lagrangian field 
shows that such realization of a free-field theory is nontrivial. While at small 
Euclidean momenta, and therefore at large distance, there is a substantial 
numerical evidence for an essentially Gaussian behavior of G{x) [132], at short 
distance renormalization group considerations lead to the asymptotic behavior 



bo and 71 are the first coefficients respectively of the /9-function and of the 
anomalous dimension of the fundamental field. We recall that a free Gaussian 
Green's function behaves like log(l/a;). Then at small distance G{x) seems 
to describe the propagation of a composite object formed by two elementary 
Gaussian excitations, suggesting an interesting hadronization picture: in the 
large- A^ limit, the Lagrangian fields U, playing the role of non- interacting 
hadrons, are constituted by two confined particles, which appear free in the 
large momentum limit, due to asymptotic freedom. 

Numerical investigations by Monte Carlo simulations of lattice chiral models 
in the continuum limit show that the large-A^ limit is rapidly approached, 
which confirms that the 1/A^ expansion, were it available, would be an effective 
predictive tool in the analysis of these models. 




(7.5) 



where A is a mass scale, and 




(7.6) 
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1.2 Principal chiral models on the lattice 



In the persistent absence of an explicit solution, the large- A?" limit of two- 
dimensional chiral models has been investigated by applying analytical and 

numerical methods of lattice field theory, such as strong- coupling expansion 
and Monte Carlo simulations. In the following we describe the main results 
achieved by these studies. 

A standard lattice version of the continuum action (7.1) is obtained by intro- 
ducing a nearest- neighbor interaction, according to Eq. (2.10): 



SU(A^) and U(A^) lattice chiral models, obtained by constraining respectively 
Ux e SU(iV) and Ux G U(iV), are expected to have the same large- A?" limit 
at fixed p. In the continuum hmit /3 — > oo, SU(A^) and U(A^) lattice actions 
should describe the same theory even at finite A^, since the additional U(l) 
degrees of freedom of U(A^) models should decouple. In other words, the U(A^) 
lattice theory represents a regularization of the SU(A^) x SU(iV) chiral field 
theory when restricting ourselves to its SU(iV) degrees of freedom, i.e. when 
considering Green's functions of the field 



e.g., 



G{x) = ^{TrUoUD, (7.9) 



whose large-distance behavior allows to define the fundamental mass M. 

At finite N, while SU(A^) lattice models should not have any singularity at 

finite j3, U(iV) lattice models should undergo a phase transition, driven by 
the U(l) degrees of freedom corresponding to the determinant of U{x). The 
determinant two-point function 

Gdix) = (det[[/t(x)C/(0)])V^ (7.10) 



behaves like x ^^1^'^^ at large x in the weak-coupling region, with f{(3,N) 
0{1/N), but drops off exponentially in strong-coupling region, where Gd{x) 
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1 , AT! 



m, = - log/? + - log ^ + (7.11) 



This would indicate the existence of a phase transition at a finite f3d in U(A^) 
lattice models. Such a transition, being driven by U(l) degrees of freedom, 
should be of the Kosterlitz-Thouless type: the mass propagating in the de- 
terminant channel should vanish at the critical point j3d and stay zero for 
larger /5. Hence for /3 > Pa this U(l) sector of the theory would decouple from 
the SU(iV) degrees of freedom, which alone determine the continuum limit 
{(3 oo) of principal chiral models. 

The large- A?" limit of principal chiral models has been investigated by Monte 
Carlo simulations of SU(A'") and U(A^) models for several large values of N, 
studying their approach to the N — oo limit [132,129]. 

Many large- A?" strong-coupling calculations have been performed which allow 
a direct study of the N — oo limit. Within the nearest- neighbor formulation 
(7.7), the large- iV strong-coupling expansion of the free energy has been cal- 
culated up to 18th order, and that of the fundamental Green's function G{x) 
(defined in Eq. (7.3)) up to 15th order [75,130]. Large- strong-coupling cal- 
culations have been performed also on the honeycomb lattice, within the cor- 
responding nearest-neighbor formulation, which is expected to belong to the 
same class of universality with respect to the critical point P = oo. On the 
honeycomb lattice the free energy has been computed up to O (/9^^), and G{x) 
upto 0(/?2°) [130]. 

Monte Carlo simulations show that SU(A'") and U(A'") lattice chiral models 
have a peak in the specific heat 

C=^^ (7.12) 
A^dT ^ ^ 

which becomes sharper and sharper with increasing A^, suggesting the presence 
of a critical phenomenon for A^ = cxd at a finite (3^- In U(^) models the peak 
of C is observed in the region where the determinant degrees of freedom are 
massive, i.e., for jS < jSd (this feature characterizes also two-dimensional XY 
lattice models [148]). An estimate of the critical couphng Pc has been obtained 
by extrapolating the position (3peak{N) of the peak of the specific heat (at 
infinite volume) to A^ ^ oo using a finite- A^ scaling Ansatz [129] 

/3peaic(A^)^/3c + cA^-^ (7.13) 



mimicking a finite-size scaling relationship. The above Ansatz arises from the 
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idea that the parameter N may play a role quite analogous to the volume in 
the ordinary systems close to the criticality. This idea was already exploited 
in the study of one-matrix models [21,48,49], where the double scaling limit 
turns out to be very similar to finite-size scaling in a two-dimensional criti- 
cal phenomenon. The finite-A^ scaling Ansatz (7.13) has been verified in the 
similar context of the large- TV Gross- Witten phase transition, as mentioned 
in Subs. 4.7. Since e is supposed to be a critical exponent associated with 
the N = oo phase transition, it should be the same in the U(A'') and SU(A^) 
models. 

The available U(iV) and SU(A^) Monte Carlo data (at = 9, 15, 21 for U(iV) 
and N = 9,15,21,30 for SU(A^)) fit very well the Ansatz (7.13), and their 
extrapolation leads to the estimates Pc = 0.3057(3) and e = 1.5(1). The inter- 
pretation of the exponent e in this context is still an open problem. It is worth 
noticing that the value of the correlation length describing the propagation in 
the fundamental channel is finite at the phase transition: ^^'^^ ~ 2.8. 

The existence of this large- A^ phase transition is confirmed by an analysis of 
the N — oo 18th-order strong-couphng series of the free energy 



F = 2/3^ + 2/3^ + 4/3^ + 19/3^ + 96/3^° + 604/3^^ 

+4036/3^^ + ^/3^^ + ^/3^« + O {p'^) , (7.14) 

which shows a second-order critical behavior: 

C-lP'^-lP-Pcr^ (7-15) 



with /3c = 0.3060(4) and a — 0.27(3), in agreement with the extrapolation of 
Monte Carlo data. The above estimates of /3c and a are shghtly different from 
those given in Ref. [147]; they are obtained by a more refined analysis based on 
integral approximant techniques [149-151] and by the so-called critical point 
renormalization method [152]. 

Green and Samuel argued that the large- phase transition of principal chiral 
models on the lattice is nothing but the large-A^ limit of the determinant 
phase transition present in U(A^) lattice models [52,111]. According to this 
conjecture, Pd and Ppeak should both converge to /3c in the large- A^ limit, and 
the order of the determinant phase transition would change from the infinite 
order of the Kosterlitz-Thouless mechanism to a second order with divergent 
specific heat. The available Monte Carlo data of U(A^) lattice models at large 
A^ provide only a partial confirmation of this scenario; one can just get a 
hint that Pd{N) is also approaching /3c with increasing A^. The large- A^ phase 
transition of the SU(A^) models could then be explained by the fact that the 
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large-A^ limit of the SU(A^) theory is the same as the large- A^" limit of the 
\J{N) theory. 

The large- character expansion of the mass propagating in the deter- 
minant channel has been calculated up to 6th order in the strong-coupling 
region, indicating a critical point (determined by the zero of the rrid series) 
slightly larger than our determination of /Sc'- l3d{N=oo) ~ 0.324 [52]. This dis- 
crepancy might be explained either by the shortness of the available character 
expansion of or by the fact that such a determination of /3c relies on the 
absence of singular points before the strong-coupling series of vanishes, and 
therefore a non-analyticity at (3c — 0.306 would invalidate all strong-coupling 
predictions for j3 > Pc- 

It is worth mentioning another feature of this large- A?" critical behavior which 
emerges from a numerical analysis of the phase distribution of the eigenvalues 
of the link operator 

L^U.UU^: (7.16) 



the N — oo phase transition should be related to the compactification of the 
eigenvalues of L [129] , like the Gross- Witten phase transition. 

The existence of such a phase transition does not represent an obstruction to 
the use of strong-coupling expansion for the investigation of the continuum 
hmit. Indeed large- Monte Carlo data show scaling and asymptotic scahng 
(in the energy scheme) even for f3 smaller then the peak of the specific heat, 
suggesting an effective decoupling of the modes responsible for the large- A^ 
phase transition from those determining the physical continuum limit. This 
fact opens the road to tests of scaling and asymptotic scahng at A^ = oo 
based only on strong-coupling computations, given that the strong-coupling 
expansion should converge for j3 < (3c- (The strong-coupling analysis does not 
show evidence of singularities in the complex /3-plane closer to the origin than 
/9c-) 

In the continuum limit the dimensionless renormalization-group invariant func- 
tion 



turns into a function A{y) of the ratio y = /Mq only, where Mq = 
and is the second moment correlation length 
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A{y) can be expanded in powers of y around y — 0: 



oo 



(7.19) 



1=2 



and the coefficients q parameterize the difference from a generahzed Gaussian 
propagator. The zero yo of A{y) closest to the origin is related to the ratio 
M'^/M^, where M is the fundamental mass; indeed yo = -M'^/M^. M'^/M^ 
is in general different from one; it is one in Gaussian models (i.e. when A{y) — 



Numerical simulations at large A^, which allow an investigation of the region 
y > 0, have shown that the large- limit of the function A{y) is approached 
rapidly and that its behavior is essentially Gaussian for y ^ 1, indicating that 
Cj <C 1 in Eq. (7.19) [142]. Important logarithmic corrections to the Gaussian 
behavior must eventually appear at sufficiently large momenta, as predicted by 
simple weak-coupling calculations supplemented by a renormalization group 
resummation: 



for p'^/Mq :$> 1 and in the large- A?" hmit. 

The approximate Gaussian behavior at small momentum is also conffimed 
by the direct estimate of the ratio M'^/Mq obtained by extrapolating Monte 
Carlo data to iV = oo. The large-iV limit of the ratio A/P/Mq is rapidly 
approached, already at = 6 within few per mille, leading to the estimate 
M^/Mq — 0.982(2), which is very close to one [132]. Large- A?" strong-couphng 
computations of M'^/Mq provide a quite stable curve for a large region of 
values of the correlation length, which agrees (within about one per cent) 
with the continuum large- A^ value extrapolated by Monte Carlo data [147]. 

Monte Carlo simulations at large values of A^ (A^ > 6) also show that asymp- 
totic scaling predictions applied to the fundamental mass are verified within 
a few per cent at relatively small values of the correlation length (^ ^ 2) and 
even before the peak of the specific heat in the so-called "energy scheme" [153]; 
the energy scheme is obtained by replacing T with a new temperature variable 
Te oc E, where E is the internal energy density. At A'^ = oo a test of asymp- 
totic scaling may be performed by using the large-A^ strong-coupling series 
of the fundamental mass. The two-loop renormalization group and a Bethe 
Ansatz evaluation of the mass/A-parameter ratio [144] lead to the following 
large- A^ asymptotic scaling prediction in the (3e scheme: 



1 + y)- 



G{P) 




(7.20) 
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V e " V4y -^'2'*^^^)' 

Pe = ^. (7.21) 

Strong-coupling calculations, where the new coupling /3e is extracted from 
the strong-coupling series of E, show asymptotic scaling within about 5% in 
a relatively large region of values of the correlation length (1.5 < ^ < 3) 
[146,147]. 

The good behavior of the large- /9-function in the f^E scheme, and therefore 
the fact that physical quantities appear to be smooth functions of the energy, 
together with the critical behavior (7.15), can be explained by the existence 
of a non- analytical zero at (3c of the /3-function in the standard scheme: 

/?i(r)=a— (7.22) 



around f3c, where a is the critical exponent of the specific heat. This is also 
confirmed by an analysis of the strong-coupling series of the magnetic suscep- 
tibility X ^G-i which supports the relations 

l^.ll^H.-..!-" (-3) 



in the neighborhood of /3c, which are consequences of Eq. (7.22) [147]. 

We finally mention that similar results have been obtained for two-dimensional 
chiral models on the honeycomb lattice by a large- strong-coupling analysis. 
In fact an analysis of the 26th-order strong-coupling series of the free energy 
indicates the presence of a large- iV phase transition, with specific heat ex- 
ponent a = 0.17, not far from that found on the square lattice (we have no 
reasons to expect that the large- A^ phase transition on the square and honey- 
comb lattices are in the same universality class). Furthermore the mass-gap 
extracted from the 20th-order strong-coupling expansion of G{x) allows to 
check the corresponding asymptotic scaling predictions in the energy scheme 
within about 10% [147]. 



7.3 The large-N limit ofSl]{N) lattice gauge theories 



An overview of the large- A^ limit of the continuum formulation of QCD has 
been already presented in Sect. 2. In the following we report some results 
concerning the lattice approach. 
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Gauge models on the lattice have been mostly studied in their Wilson formu- 
lation 



Sy^^Np J2 T^\u^,{x)U,{x + pi)Ul{x + u)Ul{x)+h.c. 



X,IJ,>V 



(7.24) 



In view of a large- A*" analysis one may consider both SU(A^) and U(A^) models, 
since they are expected to reproduce the same statistical theory in the limit 
A*" — >• oo (at fixed As for two-dimensional chiral models, SU(A^) and U(A^) 
models should have the same continuum limit for any finite N >2. 

The phase diagram of statistical models defined by the Wilson action has been 
investigated by standard techniques, i.e., strong-coupling expansion, mean 
field [77], and Monte Carlo simulations [154-156]. These studies show the 
presence of a first-order phase transition in SU(A^) models for N > A, and 
in U(A^) models for any finite N. A first-order phase transition is then ex- 
pected also in the large-A^ limit at a finite value of /?, which is estimated to 
be f3c ~ 0.38 by mean-field calculations and by extrapolation of Monte Carlo 
results. A review of these results can be found in Ref. [157]. Some speculations 
on the large- A^ phase diagram can be also found in Refs. [28,111]. The role of 
the determinant of Wilson loops in the phase transition of U(A^) gauge mod- 
els has been investigated in Ref. [Ill] by strong-coupling character expansion, 
and in Ref. [158] by Monte Carlo simulations. 

Large- A^ mean-field calculations suggest the persistence of a first-order phase 
transition when an adjoint-representation coupling is added to the Wilson 
action [116,120]. 

The first-order phase transition of SU(A'") lattice models ai N > 3 can prob- 
ably be avoided by choosing appropriate lattice actions closer to the renor- 
malization group trajectory of the continuum limit, as shown in Ref. [159] for 
SU(5). In U(A^) models the use of such improved actions should leave a resid- 
ual transition, due to the extra U(l) degrees of freedom which should decouple 
at large (3 in order to reproduce the physical continuum limit of SU(A^) gauge 
models. 

It is worth mentioning two studies of confinement properties at large A^, ob- 
tained essentially by strong-coupling arguments. In Ref. [160], the authors 
argue that deconfinement of heavy adjoint quarks by color screening is sup- 
pressed in the large- A^ limit. At N — oo, the adjoint string tension is expected 
to be twice the fundamental string tension, as implied by factorization. In Ref. 
[23], strong-coupling based arguments point out that Wilson loops in O(A^), 
U(A^), and Sp(A^) lattice gauge theories should have the same large- A^ limit, 
and therefore these theories should share the same confinement mechanism. 
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Such results should be taken into account when studying confinement mech- 
anisms. 

Studies based on Monte Carlo simulations for > 3 have not gone beyond 
an investigation of the phase diagram, so no results concerning the continuum 
limit of SU(A^) lattice gauge theories with N > 3 have been produced. Esti- 
mates of the mass of the lightest glueball, obtained by a variational approach 
within a Hamiltonian lattice formulation, seem to indicate a rapid convergence 
of the 1/N expansion [161]. 

An important breakthrough for the study of the large- limit of SU(A^) gauge 
theories has been the introduction of the so-called reduced models. A quite 
complete review on this subject can be found in Ref . [6] . 

Eguchi and Kawai [63] pointed out that, as a consequence of the large- A^ 
factorization, one can construct one-site theories equivalent to lattice YM in 
the limit N ^ oo. The simplest example is given by the one-site matrix model 
obtained by replacing all link variables of the standard Wilson formulation 
with four SU(A'") matrices according to the simple rule 

U^{x) ^ U^. (7.25) 



This leads to the reduced action 



Sek = at/? ^ Tr [u^UMlUl + h.c. 



(7.26) 



Reduced operators, and in particular reduced Wilson loops, can be constructed 
using the correspondence (7.25). In the large- A^ limit one can prove that ex- 
pectation values of reduced Wilson loop operators satisfy the same Schwinger- 
Dyson equations as those in the Wilson formulation. Assuming that all fea- 
tures of the N — oo theory are captured by the Schwinger- Dyson equations of 
Wilson loops, the reduced model may provide a model equivalent to the stan- 
dard Wilson theory at A^ = cxd. In the proof of this equivalence the residual 
symmetry of the reduced model 

C/^ ^ Z^U^, e Zn, (7.27) 



where Zjv is the center of the SU(A^) group, plays a crucial role. Therefore, the 
equivalence in the large- A?^ limit of the Wilson formulation and the reduced 
model (7.26) is actually valid if the symmetry (7.27) is unbroken. This is 
verified only in the strong-coupling region; indeed in the weak-coupling region 
the Zff symmetry gets spontaneously broken and therefore the equivalence 
cannot be extended to weak couphng [162]. 
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In order to avoid this unwanted phenomenon of symmetry breaking and to 
extend the equivalence to the most interesting region of the continuum hmit, 
modifications of the original Eguchi-Kawai model have been proposed [64,162,163]. 
The most promising one for numerical simulation is the so-called twisted 
Eguchi-Kawai (TEK) model [64,163]. Instead of the correspondence (7.25), 
the twisted reduction prescription consists in replacing 

U^{x) ^ T{x)U^T{x)\ (7.28) 

where 

T(x) = l[(r,r^ (7.29) 
and are traceless SU(A^) matrices obeying the 't Hooft algebra 

r,r^ = z^.r^r,; (7.30) 

Zfj^i, is an element of the center of the group Z^, 

^nv = exp(^i^n^^^ , (7.31) 

where n^J^y is an antisymmetric tensor with n^J^y = \ioY \i <v. are the ma- 
trices implementing the translations by one lattice spacing in the direction 
(here it is crucial that the fields C/^ are in the adjoint representation). The 
twisted reduction apphed to the Wilson action leads to the reduced action 

^TEK = A^/? E Tr \z^,U^V,VlUl + h.c] . (7.32) 

ix>v 

The correspondence between correlation functions of the large- pure gauge 
theory and those of the reduced twisted model is obtained as follows. Let 
^[[/^(x)] be any gauge invariant functional of the field Ui^{x), then 

{A[Ui,(x)])lN=oo, YM] = (^[T(x)C/^r(x)t])[jv=oo, TEK] (7.33) 



Once again the Schwinger- Dyson equations for the reduced Wilson loops, con- 
structed using the correspondence (7.28), are identical to the loop equations in 
the Wilson formulation when oo. The residual symmetry (7.27), which 

is again crucial in the proof of the equivalence, should not be broken in the 
weak-coupling region, and therefore the equivalence should be complete in this 
case. 



79 



One can also show that: 

(i) the reduced TEK model is equivalent to the corresponding field theory on 
a periodic box of size L — y/N [6] ; 

(ii) in the large- limit finite- A'^ corrections are 0{1/N'^), just as in the SU(A'^) 
lattice gauge theory. 

Moreover, since A^^ = L"^. finite- iV corrections can be seen as finite- volume 
corrections. Therefore in twisted reduced models the large-A?^ and thermody- 
namic limits are connected and approached simultaneously. 

Monte Carlo studies of twisted reduced models at large confirm the ex- 
istence of a first-order phase transition at A^ = oo located at j3c = 0.36(2) 
[164], which is consistent with the mean-field prediction f3c — 0.38 [157]. This 
transition is a bulk transition, and it does not spoil confinement. The few 
and relatively old existing Monte Carlo results obtained in the weak-coupling 
region (cfr. e.g. Rcfs. [164-166]) seem to support a rapid approach to the 
A^ ^ oo limit of the physical quantities, and arc relatively close to the cor- 
responding results for SU(3) obtained by performing simulations within the 
Wilson formulation. This would indicate that A^ = 3 is sufficiently large to 
consider the large- A?" limit a good approximation of the theory. 

We mention that hot twisted models can be constructed, which should be 
equivalent to QCD at finite temperature in the large- A^ limit (cfr. Ref. [6] for 
details on this subject). 
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